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1.  INTRODUCTION 

In  a  proposal  to  DARPA  dated  March  27,  1989,  the  following  objectives  were 
proposed  for  research  in  the  field  of  processing  of  energetics: 

1.  To  expand  the  source  codes  which  are  being  developed  for  SDIO 
(BMDOyiST  in  the  area  of  co-rotating  twin  screw  extrusion  of  highly  filled 
suspensions  to  processing  of  energetics  in  other  continuous  processors  including 
single  screw  extruders  and  oscillating  shaft  continuous  kneaders. 

2.  To  expand  the  source  codes  which  are  being  developed  for  SDIO 
(BMDOyiST  in  the  area  of  co-rotating  twin  screw  extrusion  to  incorporate  system 
components  like  dies.  The  simulation  is  to  be  relevant  to  energetic  suspensions 
which  are  generally  viscoplastic. 

3.  To  verify  the  numerical  findings  with  experimental  data. 

4.  To  develop  an  in-house  x-ray  based  diagnostic  technique  to  determine  the 
mixing  state  of  highly  filled  suspensions. 

As  will  be  shown  in  the  following,  the  objectives  (1-2)  were  met  in  conjunction 
with  viscoplastic  fluid  model  (Herschel-Bulkley  fluid)  and  the  simulation  for  single 
screw  extrusion  could  be  verified  with  experimental  data  (objective  3).  New  x-ray 
techniques  were  also  developed  to  analyze  the  degree  of  mixing  distributions  to 
satisfy  objective  4. 

II.  RESULTS 

2.1  Simulation  of  Single  Screw  Extrusion  Process 

A  non-isothermal  model  of  the  single  screw  extrusion  processing  of 
generalized  Newtonian  fluids  was  developed.  Various  temperature  dependent 
forms  of  generalized  Newtonian  fluid  constitutive  equation  representing  the 
viscoplastic  Herschel-Bulkley  fluid  and  its  simplifications  including  Bingham 
Plastic,  Power  law  of  Ostwald-de  Waele  and  Newtonian  fluids  are  applicable.  The 
model  includes  the  generally  ignored  transverse  convection  terms  of  the  equation 
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of  energy.  The  importance  of  keeping  the  transverse  convection  terms  in  the 
analysis  is  demonstrated  by  applying  the  model  and  comparison  of  findings  to 
available  experimental  results  involving  the  transverse  flow  temperature 
distributions  in  single  screw  extruders.  The  numerical  instabilities,  arising 
principally  from  the  convection  terms,  generally  encountered  in  high  Peclet 
number  extrusion  flows,  could  be  eliminated  by  the  use  of  the  Streamline 
Upwind/Petrov-Galerkin  formulation.  The  model  is  sufficiently  general  to 
accommodate  Navier's  wall  slip  at  the  wall  boundary  condition  commonly 
encountered  during  the  processing  of  energetic  suspensions. 

Complete  manuscripts  which  summarize  these  results  are  included  in  Appendix  I. 

1 .  A.  Lawal  and  D.  M.  Kalyon,  "A  Non-isothermal  Model  of  Single  Screw  Extrusion  Processing  of 
Viscoplastic  Materials  Subject  to  Wall  Slip,"  Society  of  Plastics  Engineers  ANTES  Technical 
Papers,  38  (1 992)  21 58-21 61 . 

2.  A.  Lawal  and  D.  M.  Kalyon,  “Incorporation  of  Wall  Slip  in  Non-isothermal  Modeling  of  Single 
Screw  Extrusion  Processing,"  Proceedings  of  the  First  International  Conference  on  Transport 
Phenomena  in  Processing,  S.  Guceri,  ed.,  Technomic  Publ.  Co.  (1992)  985-996. 

3.  Z.  Ji,  A.  Gotsis  and  D.  M.  Kalyon,  “Single  Screw  Extrusion  Processing  of  Highly  Filled 
Suspensions  Including  Wall  Slip,”  Society  of  Plastics  Engineers,  ANTES  Technical  Papers,  36 
(1990)  160-163. 

4.  A.  Lawal  and  D.  M.  Kalyon,  "A  Non-isothermal  Model  of  Single  Screw  Extrusion  of  Generalized 
Newtonian  Fluids,"  accepted  for  publication  in  Numerical  Heat  Transfer,  December  1992. 


2.2  Simulation  of  Oscillatory  Kneaders 

The  oscillatory  kneading  flow  was  simulated  in  conjunction  with  the  discontinuous 
cavity  flow.  Tools  of  dynamics  were  applied  to  this  flow  in  conjunction  with 
Generalized  Newtonian  Fluids,  which  are  viscoplastic  i.e.,  obey  Herschel-Bulkley 
constitutive  equation. 

The  two-dimensional  rectangular-cavity  flow  generated  by  either  continuous  or 
discontinuous  translation  of  one  or  two  walls  is  an  important  model  flow  in  efforts 
to  simulate  mixing  occurring  in  all  types  of  continuous  processors  including 
oscillatory  kneaders.  In  this  numerical  study,  the  role  of  the  viscoplasticity  of  the 
General  Newtonian  Fluids  in  affecting  their  mixing  behavior  in  the  discontinuous, 
isothermal  and  creeping  cavity  flow  is  investigated. 
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To  study  mixing  effectiveness,  various  tools  of  dynamics  including  path  line 
tracking,  Poincare'  sections  and  Lyapunov  exponents  were  employed  along  with 
the  numerical  visualization  of  the  deformations  of  discrete  blobs  of  passive 
tracers.  The  time  evolution  of  the  state  of  mixing  represented  by  the  fractional 
coverage  of  the  cavity  area  by  the  passive  tracer,  and  influenced  by  the 
development  of  chaotic  mixing  regions,  is  shown  to  be  very  different  for 
viscoplastic  fluids  in  comparison  to  Newtonian  and  shear  thinning  fluids.  Typical 
results  are  shown  in  Figures  1  and  2,  where  the  pathlines  followed  by  the  melt 
particles  in  regular  motion  and  the  Poincare'  section  for  conditions  which  lesd  into 
chaotic  dynamics  and  mixing  are  shown. 

2.3  Die  Flows 

Our  source  codes,  which  utilize  the  Finite  Element  Method  were  modified  to  allow 
the  simulation  of  die  flows.  The  methods  developed  allowed  the  incorporation  of 
the  rotation  of  the  screw  tip  immediately  preceding  the  flow  channel.  These  yet- 
to-be-published  capabilities  are  very  important  especially  for  energetic  materials 
i.e.,  viscoplastic  suspensions  which  exhibit  yield  stress.  In  the  enclosed  Figure  3 
the  typical  development  of  the  residence  time  distributions  (pathlines  followed  by 
melt  particles)  in  a  slit  die  are  shown. 

Depending  on  the  initial  location,  some  melt  particles  move  at  high  rates,  while 
others  take  longer  to  go  through  the  slit  die.  The  role  played  by  the  inclusion  of 
the  rotation  of  the  screws  is  indicated  in  the  helical  flow  at  the  entrance  to  the  die. 
The  source  code  prepared  under  DARPA  (ARPA)  funding  has  already  been  used 
in  applications  in  other  DoD  programs,  especially  in  the  design  of  an  on-line 
rheometer.  The  U.S.  Government  has  been  provided  with  a  non-exclusive 
license  on  this  invention. 

2.4  Development  of  X-ray  Based  Diagnostic  Techniques  to  Determine 
the  Mixing  State  in  Highly  Filled  Suspensions 

One  of  the  most  challenging  aspects  of  any  mixing  operation,  where  two 
or  more  identifiable  components  are  brought  together,  is  the  characterization  of 
the  state  of  the  mixture  i.e.,  the  degree  of  mixedness  or  the  "goodness"  of  mixing. 
In  the  non-diffusive  mixing  of  a  viscous  polymeric  binder  with  solid  components, 
the  complete  description  of  the  state  of  the  mixture  would  require  the 


FIG.  2. — Poincare  section  of  296  points. 


Pathlln^  followi^d  by  different  particles 


FIG.  3. — The  FEM-based  source  codes  developed  at  SIT  include 
the  capabilities  to  merge  modeling  of  die  flows  with  flow  in  the  twin 
screw  extruder.  Here,  the  velocity  distribution  and  the  pathlines 
followed  by  different  particles  close  to  the  entry  to  the  die  and  at  the 
die  are  shown. 
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specification  of  the  sizes,  shapes,  orientations  and  the  positions  of  the  ultimate 
particles  of  the  components. 

Direct  measurement  techniques  can  provide  detailed  information  on  the 
process  and  the  development  of  the  microstructure  of  the  energetic  composite. 
The  simplest  technique  which  can  be  utilized  to  determine  the  mechanisms  of 
mixing  in  a  batch  or  continuous  mixer  is  the  injection  of  a  dye  into  a  transparent 
fluid  and  then  to  follow  the  distribution  and  the  interfacial  area  growth  of  the  non- 
diffusive  tracer.  Such  techniques  were  used  earlier  to  analyze  the  mixing 
mechanisms  occurring  in  tangential  and  partially-full  co-rotating  twin  screw 
extruders.  Ottino  and  co-workers  have  employed  the  flow  visualization  technique 
to  investigate  the  mechanisms  of  mixing  in  various  time  dependent  flows,  which 
may  be  governed  by  chaotic  dynamics.  On  the  other  hand.  White  and  co-workers 
have  followed  the  spatial  distribution  of  color  incorporated  elastomers  to  gain 
insight  into  the  mixing  occurring  in  internal  mixers.  The  dispersive  mixing 
aspects,  in  conjunction  with  simple  shear  flows,  where  the  imposed  stresses  can 
be  monitored,  were  also  characterized  employing  the  cone  and  plate  flow  of  a 
transparent  oil  containing  carbon  black  agglomerates. 

If  the  materials  of  interest  are  opaque  or  if  transparent  barrel  sections 
cannot  be  built  i.e.,  as  typical  for  inert  suspensions  which  emulate  the  behavior  of 
energetic  suspensions  or  for  energetic  suspensions,  the  mechanisms  of  mixing 
can  be  studied  through  "post-mortem"  analysis.  In  this  technique,  generally  a 
distinguishable  tracer  is  added  into  a  mixer  in  a  step  or  pulse  fashion.  Upon  a 
certain  duration  of  mixing  the  mixture  is  systematically  removed  from  the  mixer 
and  sectioned  to  allow  the  investigation  of  distributive  and  dispersive  mixing 
aspects.  Tadmor  and  co-workers  have  used  this  technique  for  probing  the 
mechanisms  of  laminar  mixing  occurring  in  co-rotating  disk  processors  and  to 
follow  dispersive  mixing  of  rubber  and  carbon  black  formulations.  Gokbora  has 
utilized  pigmented  polyethylene  granules  to  investigate  the  state  of  mixedness 
occurring  in  single  screw  extrusion  processing.  Kalyon  and  co-workers  have 
employed  color  incorporated  thermoplastic  elastomers,  followed  by  computerized 
image  analysis,  to  investigate  the  distributive  mixing  of  thermoplastic  elastomers 
in  the  regular  flighted  and  kneading  disc  elements  of  twin  screw  extruders. 


On  the  other  hand,  the  rapid  advent  of  imaging  and  sensing  technology 
has  facilitated  the  introduction  of  various  powerful  techniques,  including  the 
magnetic  resonance  imaging  and  x-ray  based  techniques,  to  the  analyses  of 
opaque  mixtures.  Kalyon,  Sinton  and  co-workers  have  employed  magnetic 
resonance  imaging,  wide-angle  x-ray  diffraction  and  x-ray  radioscopy  to 
characterize  the  amount  of  air  entrained  into  composite  suspensions  with 
elastomeric  binders  during  extrusion  processing,  which  is  related  principally  to 
the  degree  of  fill  profile  in  the  extruder.  In  conjunction  with  our  DARPA  (ARPA) 
funding  the  energy  dispersive  x-ray  and  wide  angle  x-ray  diffraction  techniques 
were  adapted  and  developed  for  the  characterization  of  goodness  of  mixing  in 
concentrated  suspensions. 


EXPERIMENTAL  PROCEDURES  MATERIAL  PROCESSING  AND  SAMPLE 
PREPARATION 

The  concentrated  suspension  of  the  preliminary  demonstration  study 
consisted  of  a  polymeric  matrix  and  76.5  volume  percent  of  solids.  The  matrix 
was  hydroxyl  terminated  polybutadiene  (HTPB)  i.e.,  a  commonly  used  binder  in 
energetics.  The  fillers  were  aluminum  (Ai)  and  ammonium  sulfate  (AS)  in  powder 
form  (commonly  used  simulants  for  energetics).  Ammonium  sulfate  was  bimodal 
in  particle  size  with  seventy  five  percent  coarse  particles  and  twenty-five  percent 
fine  particles.  The  mean  particle  sizes  of  the  coarse  and  fine  fractions  of 
ammonium  sulfate  were  200  and  20  microns,  respectively.  An  APV  MP50  twin 
screw  extruder  with  a  screw  length-to-diameter  ratio  of  15  was  employed  for 
processing  of  the  suspension  samples. 

CHARACTERIZATION  OF  MIXING  INDICES 

Various  mixing  indices  which  characterize  the  goodness  of  mixing  of  the 
suspension  samples  were  determined  on  the  basis  of  the  relative  volume-fraction 
measurements  of  the  two  filler  materials  by  three  principal  analytical  methods: 
scanning  electron  microscopy  (SEM),  electron-probe  energy-dispersive  x-ray 
analysis  (EDX)  and  wide-angle  x-ray  diffractometry  (XRD).  Mixing  indices  were 
determined  at  different  scales  of  examination  by  varying  the  sample  area  through 
alteration  of  the  size  of  the  electron  and  x-ray  probes  by  several  orders  of 
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magnitude.  The  ratio  of  the  volume  fractions  of  the  two  fillers,  aluminum  (Al)  and 
ammonium  sulfate  (AS),  (t)Ai/<t>AS  .  was  utilized  as  a  measure  of  degree 
mixedness  of  the  sample. 

If  one  makes  N  measurements  of  concentration  Cj  of  one  of  the 
components,  then  the  mean  concentration  is 

N 

c-  =  l  X  C,  (1) 

/  =  1 


where  c  should  not  differ  significantly  from  <]),  the  overall  concentration  of  the 
component,  unless  a  faulty  sampling  technique  is  used.  The  difference  between 
c  and  <t)  decreases  as  the  finite  number,  N,  of  the  characterized  samples  is 
increased.  The  measured  concentration  values  of  the  minor  component  also 
depend  on  the  sample  size.  These  values  approach  the  overall  concentration  of 
the  minor  component  (|)  as  the  sample  size  is  increased. 


A  basic  measure  of  the  homogeneity  of  a  mixture  is  the  extent  to  which  the 
concentration  in  various  regions  of  the  volume  of  the  mixture  differ  from  the  mean 
concentration.  The  variance,  s^,  arising  from  the  individual  concentration,  Cj, 
measurements,  provides  such  an  index  to  quantitatively  assess  the  degree  of 
mixedness.  The  variance  is  given  by 


s 


2 


1 

{N-^) 


N 


I  {c-cf 
i  =  ^ 


(2) 


A  small  variance  value  implies  that  most  of  the  samples  yield  concentration,  Cj 
values  which  are  close  to  the  mean  c  of  all  samples,  thus  suggesting  a 
homogeneous  system.  The  maximum  variance  occurs  if  the  components  are 
completely  segregated.  Maximum  variance  is  given  by 


=  c(1-c) 


(3) 


If  the  variance  is  normalized  to  its  maximum  value  the  resulting  parameter  is 
called  the  intensity  of  segregation,  Iseg-  This  is  given  by 


I  ^  (4) 

c(1-c) 

Intensity  of  segregation,  Iseg>  is  another  important  index  in  evaluating  goodness 
of  mixing  or  distribution.  The  intensity  of  segregation  ranges  from  unity,  for 
completely  segregated  system,  to  zero,  for  a  homogeneous  system. 

In  the  demonstration,  the  mean,  variance  and  intensity  of  segregation  of 
the  relative  concentration  fractions  of  the  two  fillers,  aluminum  and  ammonium 
sulfate,  <t)Al/4>AS.  were  determined  to  assess  the  goodness  of  mixing.  The  overall 
value  of  the  relative  volume  concentration  fraction  of  aluminum  and  ammonium 
sulfate,  (t)Ai/(t)AS.  was  0.20  ±  0.02,  as  specified  during  processing. 

ANALYTICAL  TECHNIQUES 

Electron  microscopy  coupled  with  energy-dispersive  x-ray  analysis  is  a 
powerful  tool  in  the  characterization  of  microstructure  and  the  corresponding 
elemental  distribution  in  the  multi-phase  materialss.  Here  the  two  techniques 
were  utilized  in  synergy  to  obtain  both  qualitative  and  quantitative  information  on 
the  state  of  mixedness.  First  the  spatial  distribution  of  the  filler  particles  was 
determined  by  secondary-electron  microscopic  analysis.  In  order  to  distinguish 
the  aluminum  and  ammonium  sulfate  particles  from  each  other,  elemental 
mappings  of  aluminum  and  sulfur  were  carried  out  with  energy-dispersive  x-ray 
analysis  at  the  same  locations,  where  the  secondary-electron  images  were  taken. 
The  identity  of  each  particle  in  the  SEM  images  was  thus  determined  for  further 
evaluation. 

The  energy-dispersive  x-ray  analyzer  employs  a  solid  state  detector  to 
acquire  the  total  spectrum  of  characteristic  x  rays  from  0.75  to  20  keV  (or  more). 
Large  number  of  x-ray  photons  can  be  counted  at  all  energy  levels  with  the  aid  of 
a  multi-channel  analyzer.  This  allows  for  the  rapid  evaluation  of  numerous 
elements  in  the  sample.  The  intensity  of  the  characteristic  x-radiation  is 
proportional  to  the  concentration  of  the  particular  element  in  the  sampling 


volume.  Intensity  values,  therefore,  can  be  utilized  for  quantitative  analysis.  In 
this  study  the  relative  atomic  fractions  of  the  two  elements  aluminum  and  sulfur 
were  measured  by  energy-dispersive  x-ray  analysis.  Variation  of  the  relative 
atomic  fractions  of  the  two  elements  was  utilized  to  assess  the  degree  of 
distributive  mixedness.  Standard  samples  were  prepared  by  mixing  controlled 
amounts  of  aluminum  and  ammonium  sulfate  powders  in  order  to  calibrate  the 
measured  values.  According  to  this  calibration  the  relative  atomic  fraction 
proportionality  factor  kp.at.  was  1 .21 ,  i.e., 

A/^, /A/s  =  1.21  ^^1/#  (5) 

where  and  refer  to  the  atomic  concentrations  of  aluminum  and  sulfur. 
The  relative  weight  fraction  and  relative  volume  fraction  values  for  aluminum  and 
ammonium  sulfate  components  were  then  calculated  from  the  relative  atomic 
fraction  of  aluminum  and  sulfur.  The  conversion  factor  from  relative  atomic  to 
relative  weight  fraction  was  4.898  and  from  relative  atomic  to  relative  volume 
fraction  was  7.47.  Molecular  and  crystal  structure  parameters  of  the  two 
components  were  used  in  determining  these  conversion  factors.  Finally,  the 
relative  volume  fraction  proportionality  factor  kR.voi.  was  9.04,  i.e., 

I  Ns  =9.04  (6) 

Energy-dispersive  x-ray  measurements  were  carried-out  by  using  a  series 
of  successively  smaller  probe  sizes,  at  sites  chosen  systematically,  by  utilizing  a 
grid  system.  A  JOEL  840  (40A  resolution)  scanning  electron  microscope  with 
Tractor  Northern  TN  5500  energy-dispersive  x-ray  analysis  system  was  used  in 
these  studies.  The  window  thickness  was  7.5  microns.  The  probe  sizes  in  these 
experiments  were  set  at  80  mm2,  g  mm2,  1  mm2,  g.l  mm2  and  0.01  mm2  to 
determine  the  distributive  mixing  index,  i.e.,  intensity  of  segregation,  as  a  function 
of  the  scale  of  examination. 

X-ray  diffractometry  has  been  successfully  applied  for  both  qualitative  and 
quantitative  phase  analysis  in  multi  phase  materials.  The  particular  advantage  of 
x-ray  diffraction  analysis  is  that  it  discloses  the  presence  of  a  substance,  as  that 
substance  actually  exists  in  the  sample,  and  not  in  terms  of  its  constituent 
chemical  elements.  If  the  sample  contains  more  than  one  compound  or  phase 


that  constitute  the  same  chemical  elements,  all  these  compounds  are  disclosed 
by  diffraction  analysis.  Quantitative  analysis  is  also  possible,  because  the 
intensity  of  the  diffraction  pattern  of  a  particular  phase  in  a  mixture  of  phases 
depends  on  the  concentration  of  that  phase  in  the  mixture.  The  relation  between 
integrated  intensity  Ix  and  volume  fraction  (^x  of  a  phase  is  generally  nonlinear 
since  diffracted  intensity  depends  strongly  on  the  absorption  coefficient  of  the 
mixture  pm.  which  itself  depends  on  the  concentration.  For  a  two-phase  material, 
(that  with  absorption  coefficients  p-j  and  P2  for  the  two  phases)  the  absorption 
coefficient  for  the  mixture  becomes: 

fim  =  +<|J2fi2-  (7) 

The  integrated  intensity  from  phase  1  is  then  given  by; 

II  =  Ki  4)1 /pm  (8) 

where  Ki  is  a  constant  that  depends  on  the  material  and  incident  beam  used  but 
not  on  the  concentration.  The  ratio  of  intensities  from  two  phases  1  and  2, 
however,  is  independent  of  pm  and  varies  linearly  with  concentration: 

I1/I2  =  (K1/K2)  (t)i/(t)2  (9) 

In  this  study  the  ratio  relative  volume  fractions  of  aluminum  and  ammonium 
sulfate,  (1)ai/4)as.  was  calculated  from  the  integrated  intensity  ratio,  Iai/Ias-  These 
measurements  were  carried  out  utilizing  particular  crystal-planes  of  the  two  filler 
materials.  Relatively  high  number  of  crystal-plane  reflections  were  evaluated  in 
order  to  eliminate  texture  effects  and  to  increase  accuracy.  Integrated  intensities 
of  the  first  two  reflections  (1 1 1  and  200)  of  aluminum  and  the  first  five  reflections 
(002,  011,  102,  111  and  200)  of  ammonium  sulfate  were  utilized  for  calculations. 
The  value  of  the  Kai/Kas  constant  was  taken  as  0.489  for  (111)  planes.  This 
value  was  calculated  according  to  the  I/Icorundum  values  listed  by  JCPDS.  The 
contribution  of  the  absorption  of  the  matrix  to  the  absorption  coefficient  of  the 
mixture  was  less  than  3  percent  and,  thus,  it  was  negligible. 

An  automated  GE  wide-angle  x-ray  diffractometer  was  used  for  these 
studies.  Nickel  filtered  Cu  radiation  at  30  kV  and  15  mA  was  used.  The 


scanning  speed  was  one  degree  per  minute.  A  0.2  degree  receiving  slit  was  used 
in  all  runs.  The  X-ray  probe  size  was  varied  by  using  3  degree  and  0.4  degree 
primary  beam  slits  and  a  8mm  high  window.  The  X-ray  probe  sizes  used  were 
20mm  X  8mm  and  1mm  x  8mm,  respectively,  at  a  Bragg  angle,  0,  of  20  degrees. 

TYPICAL  RESULTS  AND  DISCUSSION 

In  order  to  identify  the  aluminum  and  ammonium  sulfate  particles  and  to 
carry-out  the  quantitative  analysis  of  mixing  indices,  the  elemental  mappings  of 
aluminum  and  sulfur  were  carried-out  using  energy-dispersive  x-ray  analysis  at 
the  same  locations  of  the  scanning  electron  micrographs.  These  elemental 
mappings  are  shown  in  Figures  4b  and  4c,  respectively.  As  shown  with 
connecting  arrows  in  Figure  4,  EDX  technique  can  identify  which  particles  are 
ammonium  sulfate  and  which  are  aluminum. 

The  results  of  quantitative  energy-dispersive  x-ray  analysis  of  the  relative 
volume  fractions  of  aluminum  and  ammonium  sulfate,  (1)ai/<1)as.  are  given  in  Figure 
5.  The  statistical  analysis  of  the  results  are  listed  in  Table  I.  The  variance  of  the 
relative  volume  fraction,  (l)Ai/(t)AS.  depends  on  the  probe  size.  The  largest  probe 
size,  80  mm2,  ag  expected,  gives  rise  to  the  minimum  variation  in  the  volumetric 
ratio.  With  the  moderate  probe  sizes,  which  are  in  the  1  mm2  to  1 0  mm2  range, 
the  variance  of  the  volume  ratio  remained  the  same.  At  the  relatively  smaller 
probe  sizes  of  0.1  mm2  and  0.01  mm2,  the  deviation  of  the  relative  volume 
fraction  of  the  two  solid  components  from  the  mean  increases  considerably  with 
decreasing  probe  size.  The  marked  change  in  the  variance  of  the  relative  volume 
fraction  values  between  probe  sizes  0.1  mm2  to  1  mm2  jg  jn  good  agreement  with 
the  qualitative  microstructural  features  observed  in  scanning  electron 
microscopy.  Clustering  of  the  particles  in  300  to  600  micron  diameter  regions 
should  give  rise  to  the  observed  scatter  in  the  energy-dispersive  x-ray  analysis 
results.  Probes  smaller  than  0.01  mm2  jp  gj^g  were  not  used,  as  these  would  be 
in  the  same  size  range  as  the  particle  diameter  of  ammonium  sulfate. 

The  99%  confidence  intervals  of  the  energy-dispersive  x-ray  analysis  data 
determined  according  to  Student's-t-distribution  are  also  included  in  Table  I.  The 
mean  of  the  data  for  all  probe  sizes  is  around  0.215-0.217.  These  mean  values 
agree  with  the  input  volumetric  ratio  thus  suggesting  the  absence  of  systematic 


FIG.  4. —  Microstructure  and  the  matching  elemental  mappings  of  the  filled  suspension 
taken  by  energy-dispersive  x-ray  analyzer:  (a)  scanning  electron  micrograph  at  370X 
magnification,  (b)  matching  aluminum  mapping,  (c)  matching  sulfur  mapping. 


Relative  Volume  Fraction 


0.001  0.01  0.1  1  10  100  1000 


Scale  of  Examination,  mm^ 


FIG.  5. —  Variation  of  the  relative  volume  fraction  of  the  two 
solid  components,  aluminum  and  ammonium  sulfate,  at  different 
locations  of  the  mixture  as  a  function  of  scale  of  examination 
of  energy-dispersive  x-ray  analysis. 


TABLE  I 

STATISTICAL  PARAMETERS  OF  THE  RELATIVE  VOLUME  FRACTIONS, 
<I>Al/^AS,  for  different  EDX  PROBE  SIZES 


EDX  Probe  Size 


Relative  volume  fraction, 
^Al/^’AS,  parameters 

80  mm2 

9  mm2 

1  mm2 

0.1  mm2 

0.01  mm' 

mean 

0.216 

0.215 

0.217 

0.217 

0.216 

standard  deviation 

0.007 

0.024 

0.025 

0.034 

0.069 

99%  confidence 
interval 

±0.005 

±0.014 

±0.013 

±0.018 

±0.039 

intensity  of 
segregation 


0.0003  0.0034  0.0037  0.0068  0.0281 


errors.  On  the  other  hand,  the  standard  deviation  values  increase  from  0.007  at 
the  probe  size  of  80  mm2  to  0.07  at  a  probe  size  of  0.01  mm2,  indicating  that  the 
deviation  from  the  mean  increases  by  about  one  order  of  magnitude  with  the 
decrease  in  the  probe  size.  The  confidence  intervals  again  reflect  the  same 
observation,  with  an  increasing  range  around  the  mean  as  the  probe  size  is 
reduced.  Obviously,  if  the  probe  size  were  reduced  further,  one  would  probe  only 
completely  segregated  regions,  which  would  consist  of  either  the  aluminum  or 
ammonium  sulfate  single  particles.  Which  probe  size  (scale  of  examination) 
should  be  used  for  the  characterization  of  the  goodness  of  a  given  energetic 
mixture?  The  scale  of  examination  should  be  selected  depending  on  the 
application  and  the  relevant  ultimate  properties  of  the  mixture.  For  some 
applications  in  energetics  the  sample  area  may  be  as  large  or  larger  than  the 
maximum  probe  size  used  in  this  analysis,  i.e.,  80  mm2.  However,  for  other 
applications  it  may  be  necessary  to  achieve  a  certain  degree  of  mixedness  at  a 
smaller  scale  of  examination. 

The  results  of  the  quantitative  x-ray  diffraction  analysis  of  relative  volume 
fractions  of  aluminum  and  ammonium  sulfate,  by  employing  the  relative 
integrated  intensities,  are  given  in  Figure  6.  The  statistical  analysis  of  the  results 
are  presented  in  Table  II.  With  x-ray  diffractometry,  the  mean  of  the  relative 
volume  fraction  of  aluminum  in  the  solids,  <1)ai/<1>as.  was  determined  as  0.216.  The 
analysis  was  carried  down  to  a  probe  size  of  8  mm2.  Determination  of  the  relative 
integrated  intensity,  Iai/Ias.  values  at  smaller  probe  sizes  was  found  to  be 
complicated  by  the  preferred  crystallographic  orientation  exhibited  by  the 
ammonium  sulfate  particles.  This  preferred  orientation  was  especially  evident 
with  002,  200,  or  202  planes.  As  seen  in  Figure  4,  the  ammonium  sulfate 
particles  exhibit  relatively  flat  and  angular  surfaces,  which  can  be  preferentially 
oriented  during  processing.  Since  the  samples  are  taken  from  random  locations 
and  at  random  orientations  from  the  extruded  suspension,  a  variety  of  preferred 
orientation  characteristics  were  observed  during  x-ray  diffractometry. 

The  x-ray  diffractometry  analysis  also  indicates  that  the  various  values  of 
deviation  from  the  mean  increase  with  decreasing  scale  of  examination. 
However,  there  are  differences  between  the  wide  angle  x-ray  diffraction  and 
energy-dispersive  x-ray  analysis  data  that  emanate  from  the  applicable  depth  of 
penetration.  The  electron  beam  used  in  energy-dispersive  x-ray  analysis  is 
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FIG.  6. —  Variation  of  the  relative  volume  fraction  of  the  two 
solid  components,  aluminum  and  ammonium  sulfate,  at  different 
locations  of  the  mixture  as  a  function  of  scale  of  examination 
of  x-ray  diffractometry. 


TABLE  II 


STATISTICAL  PARAMETERS  OF  THE  RELATIVE  VOLUME 
FRACTIONS,  Oai/<E)as,  FOR  DIFFERENT  XRD  PROBE  SIZES 

XRD  Probe  Size 

Relative  volume  fraction, 
$ai/^as.  parameters 

160  mm2 

8  mm2 

mean 

0.216 

0.216 

standard  deviation 

0.012 

0.029 

99%  confidence  interval 

±0.012 

±0.017 

intensity  of  segregation 


0.0009 


0.0050 


strongly  absorbed  at  the  surface  and  the  x-ray  photons  that  make  it  to  the 
detector  come  mostly  from  the  top  few  microns  depth.  This  depth  is  smaller  than 
the  particle  sizes  of  both  solid  components  of  the  formulation.  However,  in  the 
wide-angle  x-ray  diffraction  measurements  95%  of  the  information  is  derived  from 
the  top  50  to  100  microns  of  the  suspension,  depending  on  the  filler  type  and  the 
diffraction  angle  0.  The  penetration  depth  can  be  estimated  from: 

X  =  3sine /2|im  (10) 

The  penetration  depth  obtained  with  the  x-ray  diffraction  technique  is  greater  than 
the  electron  probe  measurements  and,  thus,  x-ray  diffraction  results  are  more 
representative  of  the  bulk. 

The  intensity  of  segregation  index,  Iseg,  values  were  determined  from  the 
energy-dispersive  x-ray  and  x-ray  diffractometry  data,  as  shown  in  Tables  I  and 
II.  The  intensity  of  segregation  values  increase  with  the  decreasing  scale  of 
examination  in  both  techniques.  The  intensity  values  obtained  with  energy- 
dispersive  x-ray  analysis  and  x-ray  diffraction  are  close,  indicating  that  both 
techniques  are  capable  of  characterizing  the  degree  of  distributive  mixedness  of 
the  solid  ingredients,  thus  fulfilling  the  objective  of  our  DARPA  proposal. 

A  comprehensive  manuscript  which  summarizes  the  results  (R.  Yazici  and  D. 
M.  Kalyon,  "Degree  of  Mixing  Analyses  of  Concentrated  Suspensions  by  Electron 
Probe  and  X-Ray  Diffraction,"  accepted  for  publication  in  Rubber  Chem.  and 
Techn.,  October  1 992)  is  included  in  the  Appendix  . 
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1.  INTRODUCTION 

The  fiiifie  ecRw  cjcmder  ■  one  of  ttie  nos  widely  taed  nactBnes  to  cbe  prooeusi^  In  *<**t*<?« 

to  IB  ohvKW  aobd  caneyng,  *>**»"»|  kki  prMraneaaao  ftacDooe.  u  le  ■»— h  for  nore  ipeoAlixed 

uikt  orvoived  m  muiBg  or  oogpouattag  (M}.  The  eua  of  (be  prcMi  sudy  u  lo  gtta  some 
fuwUfnfnul  angM  out  te  dyamci  of  die  dow  dsnog  (be  namg  of  wpwMMne  to  emiden.  Then 
uc  (WO  cocnpbeatiaai  rngg  frag  Ou  flow  problem.  The  dm  a  u  tbe  eoiaaoa  of  tbe 

gt  of  diffecflbel  equeaoae  cba  dracnbe  da  problan  a  (be  compbcued  geogeoy  of  da  exaiean.  The 
leoBid  oomplKgian  a  oaoBMgd  g  tte  ttaotogaal  bebeirior  of  da  ineieneL  Sun  te  (oorphotegy  of 
(be  g^wawB  cbngg  dggg  oaoDagae  gnring.  da  rbeologicel  betaevgr  of  da  maen^  eleo  ebragg 
[3].  U  is  dae  neeegwy  m  mapkmmm  •  dbev  vaeony  ttm  depatidi  <»  da  cacg  of  mixng.  fnadaraoR. 
highly  flUed  m^rmim  boro  baa  ioaa  a  slg  a  tbe  weUi  of  csgaUary  aid  raoooiat  ftiinmaat  (6). 
The  oaiAury  ao-cl^  braagry  oondraoBS  jbg.  may  kk  be  oMd  a  da  amulaab  of  flows  ovolviog 
(beg  msKfiaU.  b  fea  da  dip  gdocad  flow  cao  hronme  da  dcrotuag  mode  of  flow  ooder  oenaa  raea 
levels  (6}.  a  (be  dcgrae  dw  da  cna  daotogicel  behevgr  of  tbe  oiaenal  becomes  obecurad. 

Slip  coodmoBs  bero  bees  laad  pi'oiiaily.  a  a  owacncd  csoroaoDoe.  (nosdy  a  aaeoipa  (o  icduoe 
medamaical  aoguiatoa  aoMd  te  oaraos  aid  in  qnrodaa  swell  probleau  a  coapacoob  with  FEM 
aod  vucoelMc  added  (7  J).  Slip  a  praosue  ab  drag  adueod  isotbenul  flows  of  Newuag  flubs  a 
fibgle  anew  eatiuifcii  wa  aadad  a  (9]  uaag  ample  1<D  aod  2-D  models  aod  aaalyticb 

efpiriiiciH  io  (be  praea  papec.  coapbcaad  2*0  problems  era  rnnnflfirri  for  wbiQb  ao  aailvncai 
Mlaioos  an  avalabk. 

1 .1  Slip  at  the  Wall 

The  Blip  behavior  of  a  aapaaoD  a  da  aolid  walla  u  ebaraeanaod  a  anus  of  a  "slip  veloary".  It  is 
a  (be  afferenoe  bstwaca  da  roloory  of  (be  solid  wall  surface  (l^o)  aid  (be  vdoarr  of  da  lever 
of  tbe  fltnd  adneas  a  da  wiU  (a).  Tbe  al^  vdoory  u  generally  fotod  u  be  proporoonal  to  tbe  wall 
sbeer  aica  (e«)  [6].  a  alae  faad  a  da  aKioatruciuial  antes  of  aotaioani  of  daob^  (7).  u, 
goxral.  tbe  slip  veloaiy  can  be  proo  by  (6]: 

«  I'o  -V  -  alt.  -t^r  (1) 


wtaae  m  is  geoenlly  cloa  a  1  aal  u  (bo  vala  of  tbe  shear  rae  below  wtaieb  oe  slip  occun. 

1.2  Constitunve  Equation 

A  vismpliaic  gefnaligri  Nesranun  aodel  writb  plasoc  vacosty  j\o  oad  yield  strea  x,  wu  used  a  ibu 
study.  The  viscoaiy  fiamoD  a  giw  a  (10]; 

_  .  1  -  e*p(-«v> 

(2) 

7 

wbee  da  vuctaity  mrthnrrs  no  depends  on  da  spetafic  energy  input,  accading  a  tbe  foUowag 
eqoaaon: 

(i-«ip(-a£,)] 

and  da  speofle  energy  apa  aairta  (11]: 

DE,  dE,  o£,  3£,  1  • 

where  p  u  da  denaty.  lyl  is  da  second  avaiiaa  of  the  strain  rae.  and  I  x  I  is  da  secoad  avanaa  of  da 


(3) 


(4) 


Tbe  parameters  iv,  and  t.  are  alao  allowed  a  change  with  (emperamre  wiA  in  Anlaaa  type  of 
relanonibip: 

r{a{T)=r)o{To)tjp{-B'7<T-To}]  (5) 

%(r)  =  T,{ro)exp(-£"rfr-ro)] 


2.  THE  FINITE  ELEMENT  MODEL  OF  SINGLE  SCREW  EXTRUSION 

The  traeama  and  da  aasnmpoono  of  da  classical  tbeocy  of  sagle  screw  (l«4]  wu  ^  m  da 

preoem  sock  In  tins  daory.  da  cbaonel-Ute  dnmnns  aai  are  fomed  beiwacn  da  icrew  and  da  bairel 
may  be  unwoiad  a  fom  a  angle  atnigbi  cbranel  (U).  Then,  tbe  bairel  may  be  aasaned  a  move  wub 
respaa  a  da  screw  wnb  w  axial  snd  a  craa<<bSBDel  veloory  giyup^vvi,  eampNaslat  a  da 
angle  of  (be  screw  fli|tx.  When  tbeee  aisunpiioas  are  rnnk.  a  ctxaanmro  oqamen  of  da 
Newamao  Ruid  type  a  used  and  da  oattu  anna  and  da  gravity  are  aglsGad.  Funben&ore.  following 
da  coivenuonal  lubncaDon  usianpiion.  da  cbioga  of  velootas  a  da  down  ^iTcnon  a,  are 

neglocad.  upai  oompeiuon  a  iboa  a  i.  aod  y>  duecnons.  However,  tee  a  da  vtacoua  oargy 
dunpMon.  da  vucoesy  vatia  along  tte  diannrl  Thus,  at  erory  aecoon  along  da  cbnaei.  da  veloaiy 
dtsmbotana  depend  on  da  down  ebaanei  rtisance  t.  lo  all  ifaea  DcamoBa.  da  coDOjauiy  and 

da  luaar  momeaum  aquanone  are  solved. 


a 

0 

(6) 

d  ^  a  du. 

(T) 

d ,  ^ ,  a  du. 

ap 

dx  dx  <fy  ^ 

(8) 

^ ,  a ,  a% 

-- 

dr 

(9) 

when  da  iTis  veleaty  veear  and  P  ia  da  pmiat. 

The  pcnalty/GaleiiiD  Fiiaa  Hcinwe  Method  wa  ueed  a  da  pieaae  wak  to  (baoena  (6>-(9). 

Detmla  of  Ibis  actuque  were  giveo  u  ( 1 1  ].  A  penalry  factor  of  10*  tima  tbe  alatr  vucoeuy  wa  used 
to  da  calcnlmonA 


The  slip  law  g  eq.  (1)  is  a  combioaiioD  of  tbe  force  boundary  caq  vdocity  in  da 

ungcaial  diieeaan  of  tbe  boundary  surtace.  However,  da  normal  veloaty  on  da  tyrv  surface  u  m 
>^«wnal  They  are  qa*  difiaiit  to  couple  wnb  FEM  at  da  same  bounday.  To  solve  da  problon.  we 
pfcaea  bere  a  uavcml  tranifagatioo  of  da  esseimal  coodxooo  ato  a  oaanl  coaisaoo  form  by 
uaraduuDg  anetber  penalty  factor.  To  do  ra.  we  first  wna  botb  da  aUp  ee»vl«MB  (here  we  use  da  value 
m«l  in  eq.  1)  and  da  ipeofied  velocity  coodiboo  u  ; 


S  -5^ 


a 


(10) 


(II) 


wbere  7  and  7  are  da  tingamal  aod  da  normal  vecton  of  tbe  sitp  boundary,  lespecnvely.  A  penalty 
pasameary.  u  dan  introduced  aod  eq.  (11)  is  lewntien  m  analogy  to  eq.  (10)  for  da  normal  compoam 
N  of  da  nias  vector  ei  tbe  ai^  boondduy; 


(12) 

Eq.  (12)  transforms  tbe  earraial  coodidan  (eq.  11)  into  da  snia  fona  a  eq.  (10).  aod  it  tvyv'riMa 
eqgvaka  a  eq.  (11)  at  tbe  bmii  of  0  becanre  N  is  finiie  st  da  y  in  piacaa.  da  penalry 

paiamear  7.  taka  1  very  nnall  value.  This  tnnsfoimanao  impiia  no  matenal.  nor 

geancixy  of  da  boundary  of  da  flow  demau  sd  it  u  unveral  in  das  aena.  Of  coura.  my  kmd  of 
DOD'lioBv  form  of  dns  BanafannaDon  ts  alao  lypiicjble.  a  long  a  »  liait  ia  «q.  (11). 

The  preaenre  of  tbe  boaiduy  eondibani  a  da  fonn  of  eq.  (10)  and  (12)  ensa  da  ned  w  modify  da 
ryiam  global  menu  by  adding  to  da  appnpnaie  nodes  da  terms  fian  daa  equBoa  mvolviiig  da 
unknown  veloabes.  After  da  aoluaoD  vecar  u  found,  tbe  foUowgg  equaboo  will  pve  da  '‘slip 
vdoaty"  at  each  pond  at  tbe  waU: 


where  7*0  is  a  reference  a 


B  and  £*7  and  maeiial  panmeien. 


(13) 
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(15) 


Ncj<no(iieRntl  ctkalaaoo  n  mapkeammad  aim  te  vekxttr  eo  previow  enai  mcdob  ate  kaowa.  A 
nepwue  coaDauanoD  wtib  lespea  lo  Koorttane  amilar  to  Cnnk'Nicolaoo  oieibod  u  Med  wnb  FEM 
to  solve  Ute  eDBTBy  eqiuaao  wtacd  n  iMd: 


^7*  dT  dr  o  .  oT  5  ,  3T  . 

•  ‘■U.- - U,-; - 

a:  or  cry  ax  ax  try  ify 


(U1 


ou,  ou.  ou,  ouy  du.  du. 

“■'xi’5 - - •  -r— -‘twT” 

dr  cTv  dy  ax  dx  dy 

where  Deal  couiuctioo  aJoof  z  a  oegiecied.  The  tame  rectmique  u  applied  to  solve  tbe  epunoo  for  tbe 
speonc  eoergy  input.  Tbe  roaienai  was  amumed  to  have  ao  caergy  input  value  of  £,  of  lO^kJikg  on  tOe 
eouN' 


3.  RESULTS  AND  DISCUSSION 

The  mesb  u  the  caJcuiaooos  oescnbug  the  x*y  anas  aecoon  of  tbe  cfaannr,i  formed  m  a  angle  aow 
extruder  u  shown  m  Rg.  1.  The  operanng  candiaooi  are  pveo  to  Table  1.  The  x  std  z  componenu  of  tbe 
barrel  veloary  (V^  and  V„)  tha  mult  fzoni  a  40  rpo  manaoal  speed  of  ibe  screw  are  32  mm/iee  aod 
1 02  mm/sec  respecovelv  Tbe  rheoiopcal  and  dsennai  panmeten  to  tbe  preaeac  study  for  a  test  tntaml 
are  listed  lo  Table  0  aod  have  values  typical  of  a  suspenstoo  of  a  filler  witbis  a 

poiymenc  mainx  [11]  Veloary.  sres.  lempencun  and  specific  ezargy  profiks  were  caknlated  as  a 
fuocDOO  ot  X.  y  sod  z. 

Tbe  effea  of  tbe  slip  on  tbe  vdoaty  profile  u  dbown  to  Rg  2  Here  tbe  slip  veloaty  u  oomaltaed  over 
ihe  barrel  veloaty.  to  that  tbe  grapo  sbowt  tbe  pereeot  of  slip  as  that  wall.  The  amooni  of  elip  to  the  x 
direoiooiairye  1 )  u  ooi  effeoed  by  tbe  me  of  preaaunzaooa  b  the  down  dtaimel  Aiecnoo.  aa  «  ooiy 
refiects  tbe  shear  stress  x.  Id  tbe  oeotral  part  of  tbe  tbe  tbp  is 

ccxBtint  isbout  28%).  but  tt  greatly  near  the  ops.  abow  «*>•»  tbe  slip  akBg  tbe  aoew 

suriace  is  almo«  aeio  The  tbp  aloog  xicurve  2  aod  3)  is  sfiecied  fry  dP/dz  reflecoog  tbe  cbreigei  of  v 
aioog  X.  Agaio.  in  s  greater  poitNO  of  tbe  ccoiral  pan  of  tbe  chaooel.  tbe  tbp  is  more  toafanB  cbM  tbai  to 
x-direcQoo  .witb  s  smaller  (about  15%)  value,  with  mucb  smaller  value  (less  duo  10%)  for  die  slip  m  the 
screw  sunscetcurve  4  and  5)  However,  oear  tbe  two  ops.  tbe  sbp  is  decreased  toward  inside  cbKBd, 
which  IS  10  cooirasi  to  tbe  slip  to  z  *>^1  y  direcooos. 

z  comDooen  veiocirv  of  tbe  fluid  with  bedt  flow  lo  tbe  m  tbe  form  of  tstwels  se  shown  b  Rg.  3 

S!id  eSeci  is  cleirtv  sees  oesr  the  two  Dps.  losde  tbe  ctunoeL  (be  cootoun  of  tbe  slip  aod  o>siip  are 
Quabiativelv  same,  but  tbe  reduces  orwnsicaUy  **>»  Ttas.  of  coiwe.  anil  reduce  tbe  well 

soear  stress  level,  total  stress  leveL  aod  »*■*«  tbe  *»<Tm  energy  aod  tbe  teznperattne.  If  tbe  of  tbe 

SDCM  Bvanare  Md  toe  iprnfr  coargy  Bpre  arere  pknred  bare  for  bocb  alp  and  no  dtp  ereea.  die 
ctaaouis  wiU  ibow  dBt  abp  aaan  dare  reore  reafsB  ttBi  dare  vib  the  no  •tp  oonddua. 

TeaperHHR  prafilas  a  Pa  fiea  UmB  aa  shown  in  Hg-  4.  (a)  ia  da  eaaa  amb  tbs  no  wall  slip  caacv 
snd  (b)  u  tbe  sbp  caae.  The  eonrear  of  T  prar  erfu  to  300  K  oo  tbe  baaaL  da  arenad  arepenore  fiv 
tbe  berreL  The  ■paffc’—ww  ^  i^a  aafamrea  of  Pa  fired  a  wub  da  barel  as  da  borealary 

ctadiQoo  uaaad  of  da  oe  of  a  wafKfkbd  beat  tnnate  eoefiicaw,  may  ttatidui'r  aoree  cnw  a  date 
calculaboot.  Tbe  adubmc  cooduuna  a  da  soew  n  nth  «♦«*  Iqqi  nauance  »"** « 

reauii  a  Bgb  values  of  T  dare.  cspeaaOy  a  tbe  no-atap  caae.  The  values  of  da  saaperaore  wbeo  sip  u 
ocomog  ere  lower  tbao  **»"*»  oMwad  wia  no  sip  a  da  wall.  Be  profile  is  more  sbUobb. 

The  efiea  of  da  sbp  oondmai  oo  da  ptMaavaDon  ability  of  da  extruder  end  da  Pgandetac  of  da 
flow  rsa  00  Ua  reuoonal  epaad  at  Maaa  a  Rgreea  5  aod  6.  The  Aow  na  aod  da  pnaainiinnr 
wdueved  are  inatdepaodax  pammm  m  an  extruder  aad  Pare  two  fifota  sboald  bs  ainwinfid  agadar. 
As  I  sbnild  be  expwmid  da  fiow  om  a  sabaad  for  Pa  nae  nxsooBal  spaaL  wbsn  da  assenal  sips  m 
(be  walla.  Also,  tbe  prareunreM  obcaaad  aoder  Pe  sane  Aow  na  a  reduced  fry  da  sip.  An  eapaical 
epiiaoi  tba  fia  da  dau  repomd  bare  md  retaaea  Pc  flow  tea  and  da  pressunxmoe  for  sncb  m 
exnuder  opeiaiug  at  Q  ipm.  wbes  tbe  oreaial  alps  a  da  wall  wip  a  mrftnmi  a  U  gran  below: 


QmO^  ltr*  ( I  -0.12»o)0-0i)7»10^(l  +9*o)«~ 


By  Ibe  fwffk’anii  a  da  two  anna,  it  can  be  seen  tram  tba  dat  da  effea  of  the 

wall  sbp  a  BOR  pctreouDoeo  oo  da  efiecnve  loss  of  da  fi  iniiiisiini  ability  of  da  than  oo  tbe 

ihjcreaje  of  da  flow  rwe  for  a  pvoo  racanonal  speed.  In  odar  wards,  for  s  deand  piraiimiaiiai  na  tbe 
flow  rme  will  deoesse  dram sh rally  for  Pgb  values  of  the  sbp  wtan  da  tpm  remans  tbe 

same. 


4.  CONCLUSIONS 

A  coaprcbcnavc  Bsdieaaocal  mode)  of  tbe  singie  screw  mnnaiOD  pnoess  wre  developed  "wng  toe 
Robe  Eknaar  Method  eod  e  novel  tacbneqia  wu  wnedocnd  to  inoorperaa  da  ilip  boonday  oonPaoo 
Ml  ell  walls.  Rmbetmore.  da  shear  visocaty  matmal  hmoioa  was  varied  wob  da  energy  tmi 

BoP  sbp  aid  specific  energy  apot  dependeia  theological  behanar  are  tmts  of  tbe  proeeeBbg  of  Pghly 
filled  aopenBona.  ThM.  da  preaenad  acbmpas  Poold  be  espeasUy  oefiil  for  prePcoag  tbe 
pmmmm  bcbsvwr  of  sncb  aapeoBOOS  a  angle  soew  cjiuudai. 

The  flMdel  was  ppbed  to  a  tyfieal  flow  of  faighly  filled  suspenaaos  a  s  angle  screw  extnider.  It  was 
found  dal  tla  acctnDlaad  energy  *«*«*«■<  the  flow  field  to  a  gieaer  degree  a  da  naiial  atagec  of  da 
flow.  The  dangag  nanne  of  tbe  vuoosiiy  arts  fotad.  however,  to  beve  a  prefouod  efiba  on  oarelaove 
quananes  bkO  as  spscfik  energy  input  end  aapennire.  Oo  da  odar  baxl.  Pe  ability  of  da  aaanB  a 
sip  a  da  walls  afieca  eaeoavely  tbe  ofannng  cbarenmaics  of  da  exnoda.  ledncing  da  g—*— 
flow  me  for  a  given  prcaaiiie  drop  a  da  da  downsneam.  Oo  da  oPcr  band,  da  itap  ***<>^j«*«  tbe 
prereaizanan  abibty  of  da  exnuder.  It  ndnces  da  energy  tba  a  ecannalsud  on  tbe  ac  u  panes 

dsoBgb  da  *■»»>■*«»«  advenely  tbe  mixing  of  gjo  reduces  da 
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TABLE  1. 

Operating  Conditions  and  Geometry 
of  the  Extruder 

Screw  Diameter  =  50.8  mm 
Pitch  =  50.8  mm 

Screw  Rotational  Speed  =  40  ipm 
Inlet  Temperarure  =  300  K 
Barrel  Temperature  =  300  K 
Full  Channels 
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TABLE  2. 


Rheological  and  Thermal  Properties 
of  the  Test  Material 


T|«  =  14000  Paj 
Xj,  =  1100  Pa 
/I  =  50  j 

B  =  5  *  10"^  Jlkg 
Bt  =  0 


— a  •>  0.0000  m/MPa  • 

-  a  -  0.0092  m/MPa  s 

—  a  —  0.0020  m/MPa  a 
— a  “  0.4600  m/MPa  s 


Figure  4 
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Figure  6 
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Figure  5 
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INCORPORATION  OF  WALL  SLIP  IN  NON-ISOTHERMAL 
MODELING  OF  SINGLE  SCREW  EXTRUSION  PROCESSING 

by 


Adeniyi  Lawal  and  Dilhan  M.  Kalyon 


Highly  Filled  Materials  Institute 
Stevens  Institute  of  Technology 


INTRODUCTION 

The  single  screw  extruder  is  one  of  the  most  widely  used  continuous  processors  for 
conveying  of  solids,  melting,  pressurization,  blending,  reactive  extmsion  and  compounding  of 
a  wide  variety  of  materials.  TTie  experimental  and/or  simulation  based  knowledge  gained  over 
the  years  has  helped  our  understanding  of  the  single  screw  extrusion  process  (1,  2),  which  in 
turn  has  led  to  the  improvement  in  the  design  of  this  processing  machine. 

In  the  numerical  simulation  of  the  flow  dynamics  and  heat  transfer  occurring  in  the 
single  screw  extrusion  processing  of  a  variety  of  complicated  materials,  including  concentrated 
suspensions,  three  complications  are  manifest.  The  geometry  of  the  extruder  is  complex.  This 
can  be  effectively  dealt  with  by  using  the  finite  element  method  in  the  solution  of  the 
conservation  equations.  The  second  complication  is  related  to  the  rheological  behavior  of  non- 
Newtonian  materials.  Besides  exhibiting  shear  rate  and  temperature  dependent  material 
functions,  many  materials  also  exhibit  viscoplasticity,  i.e.  they  possess  a  limiting  stress  (yield 
stress)  below  which  they  experience  little  or  no  deformation  and  hence  flow  in  a  plug-like 
manner. 

The  third  complication  arises  from  the  behavior  of  some  viscoplastic  materials  at  solid 
boundaries.  Experimental  evidence  shows  conclusively  that  various  materials,  including 
suspensions,  gels  and  emulsions  exhibit  slip  at  the  solid  boundaries  during  flow  and 
deformation  (3-11).  The  implication  is  that  the  commonly  employed  no-slip  condition  of 
continuum  mechanics  is  not  applicable  for  a  wide  range  of  generally  viscoplastic  materials. 
Wall-slip  in  pressure  flows,  especially  the  capillary  flow  has  received  considerable  attention. 
Mooney's  (3)  investigation  of  slip  in  capillary  flow  dates  back  to  1931.  From  available 
experimental  data  on  viscometric  flows,  it  is  now  generally  accepted  that  the  slip  velocity  in 
duct  flows  can  be  described  by  the  Navier  slip  condition(12).  A  review  of  wall  slip  occurring 
at  solid  boundaries  is  available,  Schowalter  (13). 
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The  numerical  solution  of  flow  and  heat  transfer  of  polymer  melts  is  reviewed  by 
Nakazawa  et  al.  (14)  along  with  finite  element  solution  of  non-isothermal  cavity  flows  of 
Ostwald-de  Waele  model  fluids.  The  generally  encountered  highly  viscous  nature  of 
processed  materials  implies  significant  viscous  dissipation  effects  in  single  screw  extrusion.  In 
the  non-isothermal  analysis  of  single  screw  extrusion,  viscous  dissipation  and  heat  conduction 
in  the  direction  of  screw  depth  were  the  first  effects  to  be  considered  (15).  Subsequently,  heat 
conduction  in  the  other  transverse  direction  (16)  and  more  recently,  heat  convection  in  the 
down  channel  direction  (17)  were  included  in  the  solution  of  the  equation  of  energy.  Overall, 
most  non-isothermal  simulation  studies  of  single  screw  extrusion  have  ignored  the  transverse 
convection  terms,  possibly  due  to  the  fact  that  when  both  the  channel  direction  and  the 
transverse  convection  terms  are  retained  and  become  dominant,  the  commonly  utilized  standard 
Bubnov-Galerkin  method  as  well  as  its  equivalent  in  the  finite  difference  based  methods, 
produces  oscillations  in  the  numerical  solution.  This  numerical  difficulty  has  long  been 
recognized  in  the  finite  difference  and  finite  element  literature.  (18) 

The  standard  technique  for  dealing  with  spurious  oscillations  associated  with  both 
spatial  and  temporal  discretizations  of  first  derivatives  involves  the  addition  of  an  artificial 
diffusion  term  into  the  numerical  formulation  (18).  Artificial  diffusion  techniques  for  Galerkin 
finite  element  methods  can  be  accomplished  by  using  modified  forms  of  weighting  functions 
(Petrov-Galerkin)  typically  with  test  functions  which  are  one  degree  higher  than  basis 
functions  (19).  The  artificial  diffusion  term  can  also  be  incorporated  by  modifying  the  test 
functions  by  a  perturbation,  dependent  upon  the  velocity  field  (20,  21),  i.e.  Streamline 
Upwind/Petrov-Galerkin  (SUPG)  method.  For  a  parabolic  problem,  which  may  arise  from 
time  dependency  (or  in  the  case  of  steady  state  problems,  from  the  existence  of  a  predominant 
flow  direction),  the  structure  of  SUPG  method  is  only  preserved  if  time-space  elements  are 
used(22). 

In  this  study,  we  adopt  an  earlier  developed  method  of  incorporating  the  Navier  slip 
condition  for  flows  in  complex  geometries  (23).  This  method  has  been  successful  in  handling 
the  slip  condition  under  isothermal  conditions.  A  modified  Bingham  plastic  model  (24)  is  used 
to  characterize  the  rheological  behavior  of  viscoplastic  materials,  which  avoids  the  need  to 
track  yield  surfaces.  For  the  non-isothermal  analysis  in  single  screw  extruders,  the  SUPG 
formulation  is  further  developed  leading  to  an  easy  to  implement  algorithm  and  oscillation  free 
temperature  solution.  This  model  can  be  utilized  to  simulate  the  non-isothermal  processing 
behavior  of  materials  like  concentrated  suspensions,  which  exhibit  viscoplasticity  as  well  as 
wall  slip. 

ANALYSIS 


Flow  Equations 

The  governing  equations  of  the  flow  in  the  single  screw  extruder  are  given  by  the 
general  equations  of  conservation  of  mass  and  momentum,  which  are  fully  elliptic  and  three- 
dimensional.  But  the  solution  of  the  complete  set  of  three-dimensional  equations  is 
unwarranted  since  experimental  data  indicate  that  certain  simplifying  assumptions  (23)  can  be 
made  which  will  render  the  equations  parabolic  and  therefore  amenable  to  stepwise  integration 
in  the  z-direction  (primary  flow  direction)  from  prescribed  upstream  conditions.  Following  the 
parabolic  flow  an^ysis  and  introducing  the  following  dimensionless  variables; 


X 

T] 


RsX‘ 


y  =  Rg  y*  u  =  RgCO  u  *' 
R 

a'  =  :^  a*  Q  =  R^coQ* 


(la) 

(lb) 


P  =  P  * 


(1c) 
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the  equations  of  conservation  of  mass  and  momentum  become: 
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where  o)  is  the  angular  velocity  of  the  screw  and  the  asterisk  (*)  has  been  suppressed  in  the 
equations  for  convenience.  The  parabolic  nature  of  these  differential  equations  is  preserved 
through  Pm  the  mean  viscous  pressure  which  has  been  defined  for  the  primary  flow  direction, 
z  and  is  constant  for  any  x-y  plane.  We  assume  that  dpm/dz  is  constant  Equations  (3a-b)  can 
now  be  handled  separately  from  the^annel  direction  equation  (3c)  and  the  penalty  method  is 
used  in  approximating  the  pressure  p  .  The  procedure  adopted  in  this  study  for  the 
implementation  of  wall  slip  is  well  described  in  (23)  and  is  therefore  not  repeated  here. 

For  the  modified  Bingham  plastic  model  (24),  the  viscosity  T|  is  given  by: 


Xy(1  -  exp(-  n'  I  yI  )) 


Y 


(4) 


where  is  the  viscosity  coefficient,  Xy  and  n'  are  material  parameters  and  I Y I  is  defined  for 
our  case,  by: 


I  ^  I  ”  ^'^x,x  ^^y.y  ^^x,y  ^yx^  ^z,y 

Energy  Equation  and  Boundary  Conditions 

By  order-of-magnitude  analysis,  it  can  be  demonstrated  that  in  the  single  screw 
extruder,  the  transverse  convection  terms  are  not  negligible.  In  fact,  one  of  the  transverse 
convection  terms  is  generally  comparable  to  the  channel  duection  convection  term.  However, 
one  can  neglect  heat  conduction  along  the  channel  direction  and  assume  steady  state.With  these 
assumptions,  the  energy  equation  reduces  to: 
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dJ  k 


(6) 


dJ  dJ 


dx 


2dz  pC 


a 
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2T^ 


ax‘ 
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where  p  is  the  density,  Cp  the  specific  heat  and  k  is  the  thermal  conductivity.  The  fluid  enters 

the  screw  channel  at  a  temperature  and  the  barrel  temperature  is  designated  as  which  is 

different  from  the  entrance  temperature.  Variations  in  Tp  in  the  channel  direction  are  allowed. 
The  screw  surface  is  assumed  to  be  adiabatic.  Since  the  velocity  distribution  can  be 
determined  ,  this  parabolic  problem  can  be  solved  in  a  stepwise  fashion  along  the  channel 
direction  z.  Equation  6  was  solved  with  the  above  specified  boundary  conditions  in  (23)  using 
the  Bubnov-Galerkin  method,  and  the  removal  of  oscillations  in  the  temperature  solution  was 
attempted  by  a  simple  smoothing  procedure.  Here  we  implement  the  recently  developed  SUPG 
method  (which  has  rarely  been  implemented  in  polymer  processing  simulation)  by  introducing 
weighting  functions  of  the  form: 


W  =  W  +  w 


(7) 


where  W  is  a  continuous  function  of  space  and  w  is  a  discontinuous  penurbation  of  VV.  The 
weak  form  of  the  energy  equation  is  then  given  by: 


z+Az 


/  dT  3T  3T^  k  Wqy^ 


dD 


I 

•  e 


w 


aT 


ar 


ax  ay  az  ’  pC 


k  ^  .  IIL^ 


dr  dz  =  0  (8) 


where  the  last  integral  is  over  element  interiors.  The  shape  functions  are  cho.sen  to  satisfy  the 
boundary  conditions  on  the  barrel  and  must  therefore  vanish  there. 

The  shape  functions  for  the  temperature  are  similar  to  those  used  for  velocity  in  that 
they  are  linear  in  the  space  coordinates  but  the  weighting  function  W  has  a  quadratic  variation 
in  the  z-coordinate  to  improve  accuracy,  while  still  retaining  the  linear  form  in  x  and  y  (22). 
Since  the  thermal  problem  is  parabolic,  the  integration  in  z  in  Eq.  8  can  be  carried  out 
independent  of  the  other  two  directions  producing  a  part  with  known  values  of  T  (a  source 
term)  and  the  other  containing  values  of  T  being  sought.  At  this  point,  the  solution  procedure 
adopted  for  the  velocity  equations  can  then  be  applied. 

Several  approaches  have  been  suggested  in  the  literature  for  the  definition  of  the 
perturbation  function  for  SUPG.  We  have  followed  the  method  outlined  by  Hughes  et  al.  [25]. 
In  this  analysis,  an  element  in  the  arbitrary  domain  x,y,  and  z  can  be  mapped  into  an 
isoparametric  element  of  unit  dimension  in  each  of  the  transformed  co-ordinates  and 
An  element  mesh  parameter  can  then  be  defined  as 


h®  = 


2||u|| 

b 


(9) 


where  ||u||  is  the  Euclidean  norm  of  the  velocities  and  b  is  given  by 
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Hughes  et  al.  (25)  have  suggested  that  the  perturbation  function  be  taken  as: 
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■  2||u||  ^ 


VW 
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where, 

1  pCJ|u||h® 

p®  =  cothy®-  —  and  Y®  =  — 2k -  (12a,  b) 

A  value  of  o  =  1  was  used  by  Hughes  et  al.  (25).  However,  our  numerical 
experimentation  revealed  that  the  value  of  o  is  dependent  on  such  factors  as  the  element 
skewness  and  the  flow  and  heat  transfer  parameters.  Here  it  is  recommended  that  a  be 
adjusted  until  the  solution  is  oscillation  free.  For  the  range  of  parameters  considered  in  this 
study,  the  values  of  o  were  between  0.3  and  0.4  and  furnished  satisfactory  results. 


RESULTS  AND  DISCUSSION 

For  the  finite  element  solution  of  the  conservation  equations,  the  domain  is  discretized 
by  bilinear  elements  and  the  discretized  equations  for  ,  Uy  ,  ,  and  T  in  matrix  form 

were  solved  by  standard  techniques.  The  results  to  be  presented  (all  in  dimensionless  form 
defined  in  Equations  la-lc)  were  obtained  based  on  the  operating  and  processing  conditions 
listed  in  Table  1.  The  operating  conditions  are  in  the  range,  which  we  typically  utilize  in  our 
experiments  on  industrial  scale  extruders. 

Table  1. 


OPERATING  AND  PROCESSING  CONDITIONS 

q  =  1.0 

Radius  of  barrel  (R^)  = 
2.54cm 

ty  =  0.05 

Helix  angle  0  =  17.61° 

n’  =  80 

p  =  1 400  kg/m^ 

a  =  1 .0 

Cp  =  2000  J/kgK 

Screw  speed  =  40  rpm 

k  =  0.27  W/mK 

The  effect  of  wall  slip  on  the  channel' direction  velocity  U2  is  illustrated  in  Figs,  la-b 
for  a  positive  pressure  gradient  (dp^/dz  =  2.0).  The  velocity  distribution  varies  uniformly 

from  the  screw  surface  to  the  barrel  for  both  conditions  of  wail  stick  (Fig.  la)  and  wall  slip. 
Wall  slip  greatly  modulates  the  magnitude  of  the  velocity  with  the  maximum  decreasing  from 
0.953  to  0.2097  upon  wall  slip.  This  is  also  true  for  oth  values  of  positive  pressure  gradient 
not  reported  here.  As  the  value  of  dpm/dz  increase.'  i  region  of  backflow  appears,  the 
strength  and  size  of  which  increase  with  increased  slip  j;  ameter.  On  the  other  hand,  the  wall 
slip  enhances  channel  direction  velocity  for  negative  '  'ssure  gradient.  The  corresponding 
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fo  approximately  the  center  of  the  screw  surface,  thus  ensunng  that  the  fluid  feels  more 
uniformly  the  effect  of  the  barrel  temperature. 


Figure  la;  Channel  direction  velocity  distribution  for  wall  stick  condition 
(a  =  0  and  dpm/dz  =  2.0) 


Figure  lb;  Channel  direction  velocity  distribution  for  wall  slip  condition 
(a  =  1.0  and  dpm/dz  =  2.0) 


Figure  2a;  Dimensionless  temperature  distribution  at  channel  location  i  =  0.5 
for  wall  stick  condition 
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Figure  2b:  Dimensionless  temperature  distribution  at  channel  location  z  =  0.5 
for  wail  slip  condition  (a  =  1.0) 

Fig.  3a  displays  the  effect  of  wall  slip  parameter  and  pressure  gradient  on  the  xy-cross 
sectional  wall  slip  velocity  distribution.  The  region  AB  represents  the  barrel  surface  while  BA 
is  the  screw  surface.  For  the  slip  coefficient  considered,  i.e.  0.01  to  1  and  for  all  values  of 
dimensionless  pressure  gradient  dpm/dz,  the  slip  velocity  on  the  screw  surface  is  effectively 
zero. 


Figure  3a:  Slip  velocity  distribution  along  the  screw  surface  for  different 
values  of  slip  parameter  and  pressure  drop 
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The  maxima  which  occur  at  the  two  tips  of  the  screw  are  due  to  the  small  gaps  at  those 
locations.  On  the  barrel  surface,  the  slip  velocity  increases  as  the  slip  p^ameter  increases  but 
exhibits  no  appreciable  dependence  on  the  pressure  gradient,  either  positive  or  negative.  On  the 
other  hand,  the  channel  direction  slip  velocity  is  affected  by  both  the  slip  parameter  and  the 
pressure  gradient  with  a  reversal  in  behavior  occurring  from  a  change  from  negative  to  positive 
pressure  gradient  (Fig  3b). 


Figure  3b:  Channel  direction  slip  velocity  along  the  screw  surface  for  different 
values  of  slip  parameter  and  pressure  drop 

Finally,  in  Fig.4,  the  effect  of  slip  on  the  overall  pressure  drop  versus  the  volume  flow 
rate  relationship  is  depicted  for  positive  pressure  gradient  values.  Generally,  for  the  same 
pressure  gradient,  the  flow  rate  is  reduced  by  the  occurrence  of  slip  and  the  reverse  is  obtained 
for  negative  pressure  gradients. 
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Figure  4:  Pressure  drop  flow  rate  relationship  for  different  values  of  slip 
parameter 


CONCLUSION 

A  non-isothermal  model  of  the  single  screw  extrusion  processing  with  wall  slip  has 
been  developed  using  the  Streamline  Upwind/Petrov-Galerkin  (SUPG)  method.  Wall  slip 
reduces  the  capacity  of  pressurization  and  also  tends  to  produce  more  uniform  temperature 
distributions  in  the  single  screw  channel.  This  model  should  generate  a  better  predictive 
understanding  of  the  single  screw  processing  of  viscoplastic  materials  subject  to  wall  slip. 
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INTRODUCTION 

The  single  screw  extruder  Is  one  of  the  most  widely  used 
continuous  processors  for  conveying  of  solids,  melting,  pressurization, 
blending,  reactive  extrusion  and  compounding  of  a  wide  variety  of 
materials.  The  experimental  and/or  simulation  based  knowledge  gained 
over  the  years  has  helped  our  understanding  of  the  single  screw 
extrusion  process  (1,  2),  which  in  turn  has  led  to  the  Improvement  in  the 
design  of  this  processing  machine. 

In  the  numerical  simulation  of  the  flow  dynamics  and  heat  transfer 
occurring  In  the  single  screw  extrusion  processing  of  a  variety  of 
complicated  materials.  Including  concentrated  suspensions,  three 
complications  are  manifest.  The  geometry  of  the  extruder  Is  complex. 
This  can  be  effectively  dealt  with  by  using  the  finite  element  method  in 
the  solution  of  the  conservation  equations.  The  second  complication  is 
related  to  the  rheological  behavior  of  non-Newtonian  materials.  Besides 
exhibiting  shear  rate  and  temperature  dependent  material  functions, 
many  materials  also  exhibit  viscoplasticity,  i.e.  they  possess  a  limiting 
stress  (yield  stress)  below  which  they  experience  little  or  no  deformation 
and  hence  flow  in  a  plug-like  manner. 

The.  third  complication  arises  from  the  behavior  of  some 
viscoplastic  materials  at  solid  boundaries.  Experimental  evidence  shows 
conclusively  that  various  materials,  including  suspensions,  gels  and 
emulsions  exhibit  slip  at  the  solid  boundaries  during  flow  and 
deformation  (3-11).  The  implication  Is  that  the  commonly  employed  no- 
sllp  condition  of  continuum  mechanics  is  not  applicable  for  a  wide  range 
of  generally  viscoplastic  materials.  Wall-slip  In  pressure  flows,  especially 
the  capillary  flow  has  received  considerable  attention.  Mooney's  (3) 
Investigation  of  slip  in  capillary  flow  dates  back  to  1931.  From  available 
experimental  data  on  vtscometric  flows.  It  Is  now  generally  accepted  that 
the  slip  velocity  In  duct  flows  can  be  described  by  the  Navler  slip 
condltlon(12).  A  review  of  wall  slip  occurring  at  solid  boundaries  is 
available,  Schowalter(13). 

The  numerical  solution  of  flow  and  heat  transfer  of  polymer  melts 
is  reviewed  by  Nakazawa  et  al.  (14)  along  with  finite  element  solution  of 
non-lsothermal  cavity  flows  of  Ostwald-de  Waele  model  fluids.  The 
generally  encountered  highly  viscous  nature  of  processed  materials 
Implies  significant  viscous  dissipation  effects  In  single  screw  extrusion. 
In  the  non-lsothermal  analysis  of  single  screw  extrusion,  viscous 
dissipation  and  heat  conduction  In  the  direction  of  screw  depth  were  the 
first  effects  to  be  considered  (15).  Subsequently,  heat  conduction  in  the 
other  transverse  direction  (16)  and  more  recently,  heat  convection  in  the 
down  channel  direction  (17)  were  Included  In  the  solution  of  the  equation 
of  energy.  Overall,  most  non-lsothermal  simulation  studies  of  single 
screw  extrusion  have  Ignored  the  transverse  convection  terms,  possibly 
due  to  the  fact  that  when  both  the  channel  direction  and  the  transverse 
convection  terms  are  retained  and  become  dominant,  the  commonly 
utilized  standard  Bubnov-Galerkin  method  as  well  as  its  equivalent  In  the 
finite  difference  based  methods,  produces  oscillations  In  the  numerical 
solution.  This  numerical  difficulty  has  long  been  recognized  In  the  finite 
difference  and  finite  element  literature.  (18) 

The  standard  technique  for  dealing  with  spurious  oscillations 


associated  with  both  spatial  and  temporal  discretizations  of  first 
derivatives  Involves  the  addition  of  an  artificial  diffusion  term  into  the 
numerical  formulation  (18).  Artificial  diffusion  techniques  for  Galerkin 
finite  element  methods  can  be  accomplished  by  using  modified  forms  of 
weighting  functions  (Petrov-Galerkin)  typically  with  test  functions  which 
are  one  degree  higher  than  basis  functions  (19).  The  artificial  diffusion 
term  can  also  be  Incorporated  by  modifying  the  test  functions  by  a 
perturbation,  dependent  upon  the  velocity  field  (20.  21),  I.e.  Streamline 
Upwind/Petrov-Galerkin  (SUPG)  method.  For  a  parabolic  problem, 
which  may  arise  from  time  dependency  (or  In  the  case  of  steady  state 
problems,  from  the  existence  of  a  predominant  flow  direction),  the 
structure  of  SUPG  method  Is  only  preserved  If  time-space  elements  are 
used(22). 

In  this  study,  we  adopt  an  earlier  developed  method  of 
incorporating  the  Navler  slip  condition  lor  flows  in  complex  geometries 
(23).  This  method  has  been  successful  In  handling  the  slip  condition 
under  Isothermal  conditions.  A  modified  Bingham  plastic  model  (24)  Is 
used  to  characterize  the  rheological  behavior  of  viscoplastic  materials, 
which  avoids  thp  need  to  track  yield  surfaces.  For  the  non-lsothermal 
analysis  in  single  screw  extruders,  the  SUPG  formulation  is  further 
developed  leading  to  an  easy  to  implement  algorithm  and  oscillation  free 
temperature  solution.  This  model  can  be  utilized  to  simulate  the  non- 
isothermal  processing  behavior  of  materials  like  concentrated 
suspensions,  which  exhibit  viscoplasticity  as  well  as  wall  slip. 

ANALYSIS 

Flow  Equations 

The  governing  equations  of  the  flow  in  the  single  screw  extruder 
are  given  by  the  general  equations  of  conservation  of  mass  and 
momentum,  which  are  fully  elliptic  and  three-dimensional.  But  the 
solution  of  the  complete  set  of  three-dimensional  equations  Is 
unwarranted  since  experimental  data  Indicate  that  certain  simplifying 
assumptions  (23)  can  be  made  which  will  render  the  equations  parabolic 
and  therefore  amenable  to  stepwise  integration  In  the  z-directlon 
(primary  flow  direction)  from  prescribed  upstream  conditions.  Following 
the  parabolic  flow  analysis  and  introducing  the  following  dimensionless 
variables: 
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improve  accuracy,  while  still  retaining  the  linear  form  In  x  and  y  (22). 
(3a)  Since  the  thermal  problem  is  parabolic,  the  integration  In  z  In  Eq.  8  can 
be  carried  out  independent  of  the  other  two  directions  producing  a  part 
with  known  values  of  T  (a  source  term)  and  the  other  containing  values  of 
T  being  sought.  At  this  point,  the  solution  procedure  adopted  for  the 
velocity  equations  can  then  be  applied. 

Several  approaches  have  been  suggested  in  the  literature  for  the 
definition  of  the  perturbation  function  for  SUPG.  We  have  followed  the 
method  outlined  by  Hughes  et  al.  [25].  In  this  analysis,  an  element  in  the 
arbitrary  domain  x,y,  and  z  can  be  mapped  into  an  isoparametric  element 
(3c)  of  t"'!'  dimension  in  each  of  the  transformed  co-ordinates  x,  and  5.  An 
element  mesh  parameter  can  then  be  defined  as 


where  to  is  the  angular  velocity  of  the  screw  and  the  asterisk  (*)  has  been 
suppressed  in  the  equations  for  convenience.  The  parabolic  nature  of 
these  differential  equations  is  preserved  through  pm  the  mean  viscous 
pressure  which  has  been  defined  for  the  primary  flow  directon,  z  and  is 
constant  for  any  x-y  plane.  We  assume  that  dp^  /  dz  Is  constant  . 
Equations  (3a-b)  can  now  be  handled  separately  from  the  channel 
direction  equation  (3c)  and  the  penalty  method  is  used  In  approximating 
the  pressure  p  .  The  procedure  adopted  in  this  study  for  the 
implementation  of  wall  slip  is  well  described  In  (23)  and  is  therefore  not 
repeated  here. 

For  the  modified  Bingham  plastic  model  (24),  the  viscosity  t|  is 
given  by: 

_  r  yi  -exp(-n'lYl))] 


where  is  the  viscosity  coefficient,  and  n'  are  material  parameters 
and  I Y I  is  defined  for  our  case,  by: 

IyI"  =  2uf^  2ul^  +  (Ux.y  + +  “z,y  (5) 

Energy  Equatton  and  Boundary  Conditions 

By  order-of-magnitude  analysis.  It  can  be  demonstrated  that  in  the 
single  screw  extnjder,  the  transverse  convection  terms  are  not  negligible. 
In  fact,  one  of  the  transverse  convection  terms  Is  generally  comparable  to 
the  channel  direction  convection  term.  However,  one  can  neglect  heat 
conduction  along  the  channel  direction  and  assume  steady  state. With 
these  assumptions,  the  energy  equation  reduces  to: 
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where  p  Is  the  density,  Cp  the  specific  heat  and  k  is  the  thermal 
conductivity.  The  fluid  enters  the  screw  channel  at  a  temperature  T^^  and 
the  barrel  temperature  is  designated  as  Tp,  which  Is  different  from  the 
entrance  temperature.  Variations  in  Tq  in  the  channel  direction  are 
allowed.  The  screw  surface  is  assumed  to  be  adiabatic.  Since  the 
velocity  distribution  can  be  determined  ,  this  parabolic  problem  can  be 
solved  In  a  stepwise  fashion  along  the  channel  direction  z.  Equation  6 
was  solved  with  the  above  specified  boundary  conditions  in  (23)  using 
the  Bubnov-Galerkin  method,  and  the  removal  of  oscillations  in  the 
temperature  solution  was  attempted  by  a  simple  smoothing  procedure. 
Here  we  implement  the  recently  developed  SUPG  method  (which  has 
rarely  been  implemented  in  polymer  processing  simulation)  by 
Introducing  weighting  functions  of  the  form: 

W  =  W  +  w  (7 ) 

where  W  is  a  continuous  function  of  space  and  w  Is  a  discontinuous 
perturbation  of  W.  The  weak  form  of  the  energy  equation  is  then  given 
by: 
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where  the  last  integral  Is  over  element  interiors.  The  shape  functions  are 
chosen  to  satisfy  the  boundary  conditions  on  the  barrel  and  must 
therefore  vanish  there. 

The  shape  functions  for  the  temperature  are  similar  to  those  used 
for  velocity  In  that  they  are  linear  in  the  space  coordinates  but  the 
weighting  function  W  has  a  quadratic  variation  in  the  z-coordinate  to 
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where  Hull  is  the  Euclidean  norm  of  the  velocities  and  b  is  given  by 

=  (“xi  +  “yl/  *  (“xl  +  +  [^^f^  (10) 


Hughes  et  al.  (25)  have  suggested  that  the  perturbation  function  be 
taken  as: 
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A  value  of  o  =  1  was  used  by  Hughes  et  al.  (25).  However,  our 
numerical  experimentation  revealed  that  the  value  of  a  is  dependent  on 
such  factors  as  the  element  skewness  and  the  flow  and  heat  transfer 
parameters.  Here  it  Is  recommended  that  o  be  adjusted  until  the  solution 
is  oscillation  free.  For  the  range  of  parameters  considered  in  this  study, 
the  values  of  o  were  between  0.3  and  0.4  and  furnished  satisfactory 
results. 

RESULTS  AND  DISCUSSION 

For  the  finite  element  solution  of  the  conservation  equations,  the 
domain  is  discretized  by  bilinear  elements  and  the  discretized  equations 
for  u^  ,  Uy  ,  u^  ,  and  T  in  matrix  form  were  soived  by  standard 

techniques.  The  typical  geometry  of  the  single  screw  extruder  used  in  this 
study  is  shown  in  Fig.  1 .  The  results  to  be  presented  (all  In  dimensionless 
form  defined  in  Equations  la-lc)  were  obtained  based  on  the  operating 
and  processing  conditions  listed  in  Table  1.  The  operating  conditions  are 
in  the  range,  which  we  typically  utilize  in  our  experiments  on  industriai 
scale  extruders. 

The  effect  of  wall  slip  on  the  channel  direction  velocity  Uj  is 
illustrated  in  Fig.  2  for  a  positive  pressure  gradient  (dp^^/dz  =  2.0).  The 
velocity  distribution  varies  uniformly  from  the  screw  surface  to  the  barrel 
for  both  conditions  of  wall  stick  (Fig.  2a)  and  wall  slip.  Wall  slip  greatly 
modulates  the  magnitude  of  the  velocity  with  the  maximum  decreasing 
from  0.953  to  0.2097  upon  wall  slip.  This  is  also  true  for  other  values  of 
positive  pressure  gradient  not  reported  here.  As  the  value  of  dpm/dz 
increases,  a  region  of  backflow  appears,  the  strength  and  size  of  which 
increase  with  increased  slip  parameter.  On  the  other  hand,  the  wall  slip 
enhances  channel  direction  velocity  for  negative  pressure  gradient.  The 
corresponding  dimensionless  temperature  distributions  at  a  down 
channel  location  of  z  =  0.5  are  presented  in  Figs  3a-b.  The  effect  of  wall 
slip  is  essentially  the  shifting  of  the  region  of  lowest  temperature  to 
approximately  the  center  of  the  screw  surface,  thus  ensuring  that  the  fluid 
feels  more  uniformly  the  effect  of  the  barrel  temperature. 

Figure  4a  displays  the  effect  of  wall  slip  parameter  and  pressure 
gradient  on  the  xy-cross  sectional  wall  slip  velocity  distribution.  The 
region  AB  represents  the  barrel  surface  while  BA  is  the  screw  surface. 
For  the  slip  coefficient  considered,  I.e.  0.01  to  1  and  for  all  values  of 
dimensionless  pressure  gradient  dpm/dz,  the  slip  velocity  on  the  screw 
surface  is  effectively  zero.  The  maxima  which  occur  at  the  two  tips  of  the 
screw  are  due  to  the  small  gaps  at  those  locations.  On  the  barrel  surface, 
the  slip  velocity  Increases  as  the  slip  parameter  Increases  but  exhibits  no 
appreciable  dependence  on  the  pressure  gradient,  either  positive  or 
negative.  On  the  other  hand,  the  channel  direction  slip  velocity  Is  affected 
by  both  the  slip  parameter  and  the  pressure  gradient  with  a  reversal  In 
behavior  occurring  from  a  change  from  negative  to  positive  pressure 
gradient  (Fig  4b).  Finally,  in  Fig.5,  the  effect  of  slip  on  the  overall  pressure 
drop  versus  the  volume  How  rate  relationship  is  depicted  lor  positive 
pressure  gradient  values.  Generally,  for  the  same  pressure  gradient,  the 
flow  rate  is  reduced  by  the  occurence  of  slip  and  the  reverse  is  obtained 
for  negative  pressure  gradients. 
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CONCLUSION 

A  non-isothermal  model  of  the  single  screw  extrusion  processing 
with  wall  slip  has  been  developed  using  the  Streamline  Upwind  /Petrov- 
Galerkin  (SUPG)  method.  Wall  slip  reduces  the  capacity  of 
pressurization  and  also  tends  to  produce  more  uniform  temperature 
distributions  in  the  single  screw  channel.  This  model  should  generate  a 
better  predictive  understanding  of  the  single  screw  processing  of 
viscoplastic  materials  subject  to  wall  slip. 
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OPERATING  AND  PROCESSING  CONDITIONS 
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Figure  1;  Finite  element  mesh  for  the  computation  of 
the  velocity  ft  temperature  profiles. 


Figure  2b:  Channel  direction  velocity  distribution  lor 
wall  slip  condition  (osl.O  and  dPm/dzs  2.0) 


Figure  2a:  Channel  direction  velocity  diatrlbution  for 
wall  stick  condition  (a  =0  and  dPm/dz=2.0) 


Figure  3a:  Dimensionless  temperature  distribution  at 
channel  location  z=0.5  lor  wall  stick  condition 


Figure  3b:  Dimensionless  temperature  distribution  at 
channel  location  z=0.5  for  wall  slip  condition  (a=1.0) 
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Rgure  5:  Pressure  drop  flow  rate  relationship 
for  different  values  of  slip  parameter 
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Rgure  4b:  channel  direction  slip  velocity  along  the 
screw  surface  for  different  values  of  slip 
parameter  and  pressure  drop 
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ABSTRACT 


A  non-isothermal  model  of  the  single  screw  extrusion  processing  of 
generalized  Newtonian  fluids  is  presented.  Various  temperature  dependent 
forms  of  generalized  Newtonian  fluid  constitutive  equation  representing  the 
Herschel-Bulkley  fluid  and  its  simplifications  including  Bingham  Plastic,  Power 
law  of  Ostwald-de  Waele  and  Newtonian  fluids  are  applicable.  The  model 
includes  the  generally  ignored  transverse  convection  terms  of  the  equation  Of 
energy.  The  importance  of  keeping  the  transverse  convection  terms  in  the 
analysis  is  demonstrated  by  applying  the  model  and  comparison  of  findings  to 
experimental  results  involving  the  transverse  flow  temperature  distributions  in 
single  screw  extruders,  available  in  the  literature.  The  numerical  instabilities, 
arising  principally  from  the  convection  terms,  generally  encountered  in  high 
Peclet  number  extrusion  flows  could  be  eliminated  by  the  use  of  the  Streamline 
Upwind/Petrov-Galerkin  formulation.  The  model  is  sufficiently  general  to 
accommodate  Navier's  wall  slip  at  the  wall  boundary  condition  commonly 
encountered  during  the  processing  of  gels,  and  concentrated  suspensions. 

INTRODUCTION 

The  single  screw  extruder  is  one  of  the  most  widely  used  continuous 
processors  for  conveying  of  solids,  melting,  pressurization,  blending,  reactive 
extrusion  and  compounding  of  a  wide  variety  of  materials.  The  experimental 
and/or  simulation  based  knowledge  gained  over  the  years  has  helped  our 
understanding  of  the  single  screw  extrusion  process  and  are  summarized  in 
various  sources  (1-5). 

The  numerical  solution  of  flow  and  heat  transfer  of  polymer  melts  is 
reviewed  by  Nakazawa  et  al.  (6)  along  with  the  finite  element  solution  of  non- 
isothermal  cavity  flows  of  Ostwald-de  Waele  model  fluids.  The  generally 
encountered  highly  viscous  nature  of  processed  materials  implies  significant 
viscous  dissipation  effects  in  single  screw  extrusion.  Fenner  investigated 
various  solution  methods  available  for  the  flow  and  deformation  occurring  in  the 
single  screw  extruder  and  compared  the  findings  against  experimental  results 
(2).  The  geometry  of  the  single  screw  extruder  is  usually  unwound  in  the 
channel  direction  and  the  curvature  effects  are  neglected.  Zamodits  and 
Pearson  (3)  assumed  the  channel  geometry  to  be  rectangular  but  the  effects  of 


the  flight  on  the  velocity  distribution  were  ignored.  Screw  characteristic  curves 
of  volume  flow  rate  versus  pressure  generated  for  different  values  of  power  law 
index  were  presented.  Denson  and  Hwang,  using  the  Finite  Element  Method, 
analyzed  the  flow  in  several  actual  channel  geometries  but  the  fluid  was 
assumed  to  be  Newtonian.  An  equation  relating  the  throughput  to  the  pressure 
rise  was  developed  (5). 

In  the  non-isothermal  analysis  of  single  screw  extrusion,  viscous 
dissipation  and  heat  conduction  in  the  direction  of  screw  depth  were  the  first 
effects  to  be  considered  (7).  Subsequently,  heat  conduction  in  the  transverse 
direction  (8)  and  more  recently,  heat  convection  in  the  down  channel  direction 
(9)  were  included  in  the  solution  of  the  equation  of  energy.  Overall,  most  non- 
isothermal  simulation  studies  of  single  screw  extrusion  have  ignored  the 
transverse  convection  terms,  possibly  due  to  the  fact  that  when  both  the  channel 
direction  and  the  transverse  convection  terms  are  retained  and  become 
dominant,  the  commonly  utilized  standard  Bubnov-Galerkin  method  as  well  as 
its  equivalent  in  the  finite  difference  based  methods  (i.e.  the  central  difference 
scheme),  produces  oscillations  in  the  numerical  solution.  This  numerical 
difficulty  has  long  been  recognized  in  the  finite  difference  and  finite  element 
literature  (10). 

The  standard  technique  for  dealing  with  spurious  oscillations  associated 
with  both  spatial  and  temporal  discretizations  of  first  derivatives  involves  the 
addition  of  an  artificial  diffusion  term  in  the  numerical  formulation  (1 0).  Artificial 
diffusion  techniques  for  Galerkin  finite  element  methods  can  be  accomplished 
by  using  modified  forms  of  weighting  functions  (Petrov-Galerkin)  typically  with 
test  functions  which  are  one  degree  higher  than  basis  functions  (11).  The 
artificial  diffusion  term  can  also  be  incorporated  by  modifying  the  test  functions 
by  a  perturbation,  dependent  upon  the  velocity  field  (12,  13),  i.e.  Streamline 
Upwind/Petrov-Galerkin  (SUPG)  method.  For  a  parabolic  problem,  which  may 
arise  from  time  dependency  (or  in  the  case  of  steady  state  problems,  from  the 
existence  of  a  predominant  flow  direction),  the  structure  of  SUPG  method  is  only 
preserved  if  time-space  elements  are  used  (14).  In  the  finite  difference  literature, 
the  spurious  oscillations  are  usually  remedied  by.  using  a  combination  of 
central-difference  and  upwind-difference  formulations  for  the  first  order 
derivatives  (15),  a  procedure  akin  to  the  classical  upwinding  technique  of  the 
.  finite  element  method. 


In  the  following,  we  present  a  new  model  of  non-isothermal  extrusion  of 
generalized  Newtonian  fluids  in  single  screw  extruders.  In  the.  model,  the 
generally  ignored  transverse  convection  terms  of  the  equation  of  energy  are 
included,  in  conjunction  with  the  Streamline  Upwind/Petrov-Galerkin 
formulation  of  the  Finite  Element  Method.  The  modeling  results  are  compared 
with  the  experimental  results  available  in  the  literature  to  demonstrate  the 
importance  of  keeping  the  transverse  convection  terms.  The  model  also  has  the 
facility  to  include  the  Navier's  wall  slip  boundary  condition  commonly 
encountered  with  a  variety  of  fluids  (16*22).  Various  simplifications  of  the 
Herschel-Bulkley  viscoplastic  fluid  model,  i.e.,  the  Bingham  Plastic,  the  Power- 
law  fluid  of  Ostwald-de  Waele  and  the  Newtonian  fluid  can  be  accommodated. 


ANALYSIS 
Flow  Equations 

The  governing  equations  of  the  flow  in  the  single  screw  extruder  are 
given  by  the  general  equations  of  conservation  of  mass  and  momentum,  which 
are  fully  elliptic  and  three-dimensional.  But  the  solution  of  the  complete  set  of 
three-dimensional  equations  is  unwarranted  since  experimental  data  indicate 
that  certain  simplifying  assumptions  can  be  made  which  will  render  the 
equations  parabolic  and  therefore  amenable  to  stepwise  integration  in  the  z- 
direction  (primary  flow  direction)  from  prescribed  upstream  conditions.  These 
assumptions  are  that: 

(i)  the  inertia  effect  is  negligible  because  of  the  highly  viscous  nature  of  the 
generally  polymeric  materials  used  in  extrusion  and  the  creeping  nature  of 
the  resulting  flow, 

(ii)  the  flow  Is  steady, 

(iii)  the  screw  channel  can  be  unwound  in  a  helical  direction  of  the  screw, 

(iv)  the  velocity  derivatives  along  the  channel  are  small  compared  with  those 
in  the  transverse  directions  except  close  to  the  entrance,  and 

(v)  the  barrel  moves  relative  to  the  screw  at  an  angle  (i.e.  the  helix  angle)  to 
the  down  channel  direction  with  a  constant  velocity  given  by  the  screw 
rotational  speed  and  barrel  diameter. 


The  parabolic  flow  analysis  is  followed  and  the  following  dimensionless 
variables  are  introduced: 

x  =  RgX*,y=Rsy*,z=RsZ*.u  =  RsCou*  {1a) 

T|  =  moco'^~^'n*.  Q  =  R^cdQ*  ("•  b) 

p  =  moCo'^p*,  Pm  =  mo“"pm‘ 

where  Rs  is  the  screw  radius,  co  is  screw  rotational  speed,  and  n  are 
material  parameters  and  Q  is  the  volumetric  flow  rate.  The  equations  of 
conservation  of  mass  and  momentum  become: 


where  the  asterisk  (*)  denoting  the  dimensionless  forms  of  the  variables  has 
been  suppressed  in  the  equations  for  convenience.  The  parabolic  nature  of 
these  differential  equations  is  preserved”  through  pm,  the  mean  viscous  pressure 
which  has  been  defined  for  the  primary  flow  direction,  z,  and  is  constant  for  any 
x-y  plane.  The  down-channel  direction  pressure  gradient,  dpm/dZi  which  is  a 
function  of  the  channel  direction  coordinate,  z,  for  temperature  dependent  shear 
viscosity,  is  determined  by  the  requirement  that  the  conservation  of  mass 
constraint  on  the  total  volumetric  flow  rate  is  satisfied,  i.e.: 

Q- j  UzdD  =  0 
D 


(4) 


where  Q  is  set  by  the  specified  inlet  vaiue  of  dpm/dz.  The  pressure  p  ,  is 
ailowed  to  vary  in  the  transverse  directions  in  such  a  way  that  the  continuity 
equation  (Eq.  2)  is  satisfied.  Even  though  uz  varies  along  the  down-channel 
direction  z,  it  is  assumed  in  Eq.  2  that  step  sizes  in  z  direction  are  so  small  that 
3uz/5z  is  negligible  compared  with  the  other  two  terms.  Equations  3a-b  can  now 

be  handled  separately  from  the  channel  direction  equation  3c  and  the  penalty 
method  is  used  in  approximating  the  pressure  ~p. 

Boundary  Conditions  and  impiementation  of  Waii  Slip 

Applying  Bubnov-Galerkin's  method  to  Eqs.  3a-b  leads  to  the  following 
dimensionless  residual  equations: 

J  Vtc  •  TdD  -  J  ;c(n  •  1)  dr  =  0 
D  -  r  ~ 

s 

where  n  is  the  weighting  function  and  T ,  the  total  stress  tensor  is  as 
the  parabolic  flow  assumption.  On  the  boundaries,  the  traction  n 
decomposed  into  the  tangential  and  normal  components  to  obtain: 

Jv;t-TdD-  j  7c(nt :  T)t dPg  -  J  7t(nn  r^ndP  =  0  (6) 

Dr  r- 


(5) 

modified  by 
•  T  can  be 


where  t  and  n  are  respectively  the  unit  tangent  and  unit  outward  normal  vectors 
to  the  solid  boundary.  Eqs.  3a-b  are  required  to  be  satisfied  on  the  boundary. 
The  conditions  of  wall  slip  and  no  material  exchange  between  the  fluid  and  the 
boundary  in  dimensioniess  form,  can  be  expressed  as: 


‘•(“•“solid’  =  P*("‘=I’ 


"•‘“-“solid’ =  ° 


(7a) 

(7b) 


where  Ug  is  the  solid  boundary  velocity,  u  the  fluid  velocity,  (3*  =  (moO)^“VRs)13 
is  the  dimensionless  Navier's  slip  coefficient,  i.e.,  p  =  0  gives  no  slip  and  p  =  oo 
is  perfect  slip.  Eq.  7a  is  the  Navier's  slip  condition  as  applicable  to  the  two- 
dimensional  transverse  plane  (23). 


Referring  now  to  Eq.  6,  for  the  implementation  of  the  boundary 
conditions,  one  observes  that  while  the  slip  condition  easily  blends  with  it  by 
providing  the  tangential  component  of  the  traction,  the  normal  component 
presents  a  difficulty  as  it  is  not  recoverable  from  the  second  condition  Eq.  7b, 
which  gives  only  the  normal  velocity.  To  resolve  this,  the  essential  condition  of 
normal  velocity  is  transformed  into  a  mixed  condition  (24): 


nt  ‘T  =  ^•(^~^solid) 

P* 


(8a) 


nn  :T  =  ^solid) 

A-s 

where  a  new  penalty  parameter,  has  been  introduced.  As  0,  Eq.  8b 

becomes  equivalent  to  Eq.  7b  since  the  normal  component  of  the  traction  must 
be  finite  at  the  boundary.  In  practice,  a  Xg  value  in  the  range  10*''°  - 10*''° 

produces  a  normal  velocity  difference  of  the  same  order  of  magnitude. 
Equations  8a-b  are  now  substituted  Into  Eq.  6  to  give: 


JVtcT  dD  -  J  „i±L^S0|id)  J  ,iL(!izMd)ndrs  =  0  (9) 

D  r3  P 


For  the  channel  direction  momentum  equation,  the  normal  component  of 

traction  does  not  exist  which  makes  the  incorporation  of  wall  slip  by  the  Galerkin 
method  straightforward.  The  components  of  Ug^ij^j  on  the  screw  root  surface  are 

all  zero  while  on  the  barrel,  the  solid  boundary  velocity  is  given  in 

dimensionless  form  as: 

» 

“solid  =  'sin0  +  k  cos  0  (10) 

where  0  is  the  helix  angle  of  the  screw  and,  I  and  k  are  unit  vectors  in  x-  and  z- 
directions  respectively. 


Constitutive  Model 


8 


The  fluid  is  assumed  to  be  a  purely  viscous  fluid  which  can  be 
characterized  rheologicaliy  by  the  generalized  Newtonian  fluid  constitutive 
relation: 


l  =  -n(7)  Y 


(11) 


where  ti  is  the  shear  viscosity  material  function,  y  is  the  rate-of-deformation 

■» 

tensor,  and  y  is  the  deformation  rate.  The  dependence  of  the  shear  viscosity 
material  function  bn  the  deformation  rate  is  assumed  to  follow  the  modified 
Herschel-Bulkley  model,  which  for  a  temperature  dependent,  viscosity  is  given 
by: 


Tl  = 


molYl' 


,n-?  ^  Ty(f-exp(-nb|y|)) 


Y 


exp(-c*{T-To)) 


(12) 


where  n  Is  a  material  parameter  which  governs  the  sensitivity  of  the  fluid  to 
deformation  rate,  Xy,  the  apparent  yield  stress,  nb,  a  stress  growth  exponent,  c', 
the  temperature  coefficient  of  viscosity,  with  To,  the  entrance  temperature  as  the 
reference  temperature.  The  magnitude  of  the  rate-of-deformation  tensor  |y|  is 
defined  for  our  case,  by: 

IyP  =  2ul^  +  2u^j,  +  +  (13) 

where  comma  indicates  differentiation.  If  we  introduce  the  following 
dimensionless  variables: ' 


|y|  =  coIyI*,  Ty  =  moco'^xy ,  n^,  =  n^  /  co 


then  Eq.  12,  in  dimensionless  form,  becomes: 


Tl*  = 


,  <(1-exp(-nb|Y|-)) 


(14) 


exp(-c'(T-To)) 


(15) 


Various  simplifications  of  the  modified  Herschel-Bulkley  model  include  the 
Ostwald-de  Waele  power-law  fluid  {ty  =  0),  the  Newtonian  fluid  (n=1,  ty  =  0), 

the  modified  Bingham  plastic  (Papanastasiou  (25),  (n=1)),  and  the  classical 
Bingham  plastic  fluid  above  the  yield  point  {n=1,  =  «>). 


Energy  Equation  and  Boundary  Conditions 

By  an  order-of-magnitude  analysis  of  the  energy  equation,  it  can  be 
demonstrated  that  in  the  single  screw  extruder,  the  importance  of  the  transverse 
convection  terms  (pCpUx3T/3x,  pCpUy3T/3y  )  relative  to  heat  conduction  in  y- 

direction  scales  as  Peu^H^/W  where  Pe  =pCpCDRg/k  is  the  Peclet  number,  and 

H  and  W  are  respectively,  the  dimensionless  maximum  depth  and  width  of  the 
screw  channel.  Polymeric  liquids  generally  have  low  thermal  conductivity  and, 
moderately  high  values  of  p  and  Cp,  hence  high  Peclet  number.  Therefore,  the 

transverse  convection  term  will  be  negligible  only  for  small  H  and,  large  W  (i.e. 
low  aspect  ratio  geometries),  and  also  for  small  helix  angle.  If  the  transverse 
convection  term  is  compared  with  the  down  channel  convection  term,  the 
transverse  convection  term  will  be  negligible  if  the  ratio  uJ^z’/u^W  is  small, 
where  ujj  and  Uz  are  dimensionless.  From  the  continuity  equation,  and  the 
choice  of  step  size,  u^/z*  is  small,  hence  transverse  convection  can  be  ignored 
only  if  W  is  very  large  or  the  helix  angle  very  small,  the  same  conditions  as 
earlier  obtained.  Since  these  conditions  are  rarely  used  in  the  design  and 
operation  of  realistic  extruders,  the  transverse  convection  terms  should  be 
included  in  the  energy  equation  for  the  analysis  of  extrusion  flows. 

By  a  similar  analysis,  the  quantity  Peu^z*  governs  the  relative 
importance  of  the  down  channel  convection  to  heat  conduction  in  the  down 
channel  direction.  Therefore,  except  for  low  volume  flow  rates,  and  the  region 
very  close  to  the  entrance,  one  can  neglect  heat  conduction  in  the  down 
channel  direction  since  the  Peclet  number  is  generally  high.  With  these 
assumptions  and  the  assumption  of  steady  state  conditions,  using  the 
dimensionless  temperature  8  =  (T  -  To)  /  (Tb  -  To),  the  energy  equation  in 
dimensionless  form  becomes: 
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de  de 


a^e' 

^ax^'^ay^^ 


+  Gti 


(16) 


where  G=m„co""‘'(coRj2/k(Tu-  TJ  is  the  Griffith  number,  (generally  referred  to 

as  the  Brinkmann  number  for  the  Newtonian  fluid)  and  the  asterisk  (*)  denoting 
the  dimensionless  forms  has  been  suppressed  for  the  velocity  components  and 
co-ordinates.  The  Griffith  number  indicates  the  relative  importance  of  viscous 
dissipation  effects  to  heat  conduction. 

The  fluid  enters  the  screw  channel  at  a  temperature  T^.  The  barrel 
temperature  is  designated  as  T^.  and  is  different  from  the  entrance  temperature. 

Variations  in  Tb  along  the  channel  direction  are  allowed.  Although  the  model 
can  be  used  in  conjunction  with  any  type  of  thermal  boundary  condition,  we  will 
restrict  our  demonstration  here  to  the  use  of  constant  barrel  temperature  and 
adiabatic  screw  surface.  Thus,  in  dimensionless  form,  the  Initial  and  boundary 

conditions  become: 


at  z  =  0, 0  =  0 
on  barrel  surface  0  =  1 

on  screw  surface  90  /  9n  =  0  (1 7a-c) 

where  n  is  the  unit  vector  normal  to  screw  surface.  Since  the  velocity  distribution 
can  be  determined,  this  parabolic  problem  can  be  solved  in  a  stepwise  fashion 
along  the  channel  direction,  z.  Equation  16  was  solved  with  the  above  specified 
boundary  cpnditions  in  (24)  using  the  Bubnov-Galerkin  method,  and  the 
removal  of  oscillations  in  the  temperature  solution  was  attempted  by  a  simple 
smoothing  procedure.  Here  we  implement  the  recently  developed  SUPG 
method  (12-14),  (which  has  rarely  been  implemented  in  polymer  extrusion 
simulation)  by  introducing  weighting  functions  of  the  form: 

=  W^  +  w^ 

Where  W<  is  a  continuous  function  of  space  and  wt  is  a  discontinuous 
perturbation  of  W^.  The  weighted  residual  form  of  the  energy  equation  is  then 

given  by: 


(19) 


The  shape  functions  are  chosen  to  satisfy  the  boundary  conditions  on  the  barrel 
and  must  therefore  vanish  there. 


Several  approaches  have  been  suggested  in  the  literature  for  the 
definition  of  the  perturbation  function  for  SUPG.  We  have  followed  the  method 
outlined  by  Hughes  et  al.  (26).  In  this  analysis,  an' element  in  the  arbitrary 
domain  x,y,  and  z  can  be  mapped  into  an  isoparametrio  element  of  unit 
dimension  in  each  of  the  transformed  serendipity  coordinates  x,  and  C.  An 
element  mesh  parameter  can  then  be  defined  as 


(20) 


where  ||u||  is  the  Euclidean  norm  of  the  velocities  and  b  is  given  by 

=  (“x  i  +  “y  +  (‘'x  ^  +  [‘z  i]^  (21) 

Hughes  et  al.  (26)  have  suggested  that  the  perturbation  function  be  taken  as: 


,  ah®e®  _ _ 


,  (22) 


where, 

e®  =  cothv®-  and  7®  =  Pe|lu||h6/2, 

r 


(23a,  b) 


NUMERICAL  SOLUTION 


Flow  Equations 


For  the  finite  element  solution  of  the  momentum  equations,  the  domain  is 
discretized  by  bilinear  elements  in  x  and  y,  and  trial  functions  are  chosen  such 
that: 

Nod 

X®  =  X  x,®7cf(x,y)  Xj®,7f  €  D®  (24  a,  b) 

i  =  1 

where  D®  is  an  element  domain,  Nod  is  the  number  of  nodes  per  element,  and 
nf  is  the  trial  function  with  x®  representing  the  velocity  vector.  The  discretized 
equations  take  the  matrix  form: 

=  ^  (25) 

where  A  is  the  stiffness  matrix  and  f  is  a  force  vector  resulting  from  specified 

nodal  values  or  source  term  and  the  nodal  unknowns  are  recovered  from  the 
vector  X. 


Energy  Equation 

•v- 

The  shape  functions  for  the  temperature  are  similar  to  those  used  for 
velocity  in  that  they  are  linear  in  the  space  coordinates.  The  weighting  function 
Wf  has  a  trilinear  dependence  on  the  x-,  y-,  and  z-  coordinates.  Since  the 
thermal  problem  is  parabolic,  the  integration  in  z  in  Eq.  16  can  be  carried  out 
independent  of  the  other  two  directions  producing  a  part  with  known  values  of  0 
(source  term)  and  the  other  containing  values  of  0  being  sought.  At  this  point, 
the  solution  procedure  adopted  for  the  velocity  equations  can  then  be  applied. 


The  element  mesh  parameter  h®  ,  and  the  cell  Peclet  number  7®  in  Eqs. 
20  and  23  were  evaluated  based  on  the  velocity  at  the  centroid  of  the  element. 
A  value  of  a  =  1  in  Eq.  22  was  suggested  by  Hughes  et  al.  (26).  However,  our 
numerical  experimentation  revealed  that  the  value  of  a  is  dependent  on  such 

factors  as  the  element  skewness  and  the  flow  and  heat  transfer  parameters. 
Here  it  is  recommended  that  a  be  adjusted  until  the  solution  is  oscillation  free. 
For  the  range  of  parameters  considered  in  this  study,  the  values  of  a  were 
between  0.5  and  3.0  and  furnished  satisfactory  results. 


Solution  Procedure 

At  the  inlet  of  the  channel,  the  temperature  distribution  is  known,  hence 
only  the  x-,  y-,  and  z-  velocity  components  are  computed,  based  on  a  specified 
value  of  the  inlet  pressure  gradient  dpm'/dz*.  The  Picard  iteration  (27)  is  then 
-utilized  for  the  shear  rate  dependent  shear  viscosity.  After  convergence,  the 
volumetric  flow  rate  is  calculated.  We  proceed  down  the  channel  and  as  a  first 
step,  solve  for  the  temperature  distribution.  With  a  viscosity  distribution  based 
on  this  temperature  profile  and  the  value  of  dpmVdz*  at  the  previous  channel 
location,  for  the  second  step,  we  solve  for  the  velocity  components  and 
calculate  the  volumetric  flow  rate,  A  new  value  of  dpmVdz*  is  obtained  using  the 
secant  method  with  the  residual  taken  as  the  difference  between  the  calculated 
volumetric  flow  rate  and  that  obtained  at  the  inlet.  This  completes  the  second 
step  which  is  repeated  until  convergence  on  account  of  shear  rate  dependent 
shear  viscosity. 

Subsequently,  we  solve  for  a  new  temperature  distribution  and  go  back 
to  the  second  step  since  the  shear  viscosity  material  function  is  temperature 
dependent.  This  iterative  procedure  is  continued  until  the  volumetric  flow  rate 
constraint  (Eq.  4)  is  satisfied  and  simultaneously,  the  Euclidean  norms  of  the 
difference  between  successive  iterates  of  u,  and  0  are  small,  usually  ^  10-3. 
Another  step  size  is  taken  with  the  values  of  the  variables  at  the  previous 
channel  location  serving  as  the  initial  guesses.  Convergence  on  the  inner 
iteration  for  shear  rate  dependent  shear  viscosity  is  usually  achieved  in  4  to  8 
iterations  while  the  outer  Picard  iteration  on  account  of  the  temperature 
coefficient  of  viscosity  may  require. up  to  10  iterations,  with  both  values 


depending  on  the  channel  location  and  the  value  of  the  temperature  coefficient 
of  viscosity.  When  the  backflow  is  significant  but  not  severe,  the  SUPG  method 
provides  reliable  and  accurate  results  when  compared  with  the  upwind 
marching  scheme  of  finite  difference  but  the  required  number  of  iterations  goes 
up  slightly.  The  allowable  step  size  Az,  is  dictated  by  the  stability  of  the 
discretized  temperature  equation  which  is  dependent  on  the  thermal  and  flow 
parameters.  The  dimensionless  step  size  used  in  this  study  is  1.0. 


RESULTS  AND  DISCUSSION 

The  code  was  initially  tested  by  comparison  of  the  well-known  results 
available  in  the  literature  for  the  square  cavity  flow  of  a  Newtonian  fluid  in 
conjunction  with  the  no-slip  at  the  wall  boundary  condition.  As  indicated  in  Fig. 
1  a,  excellent  agreement  was  found  between  our  results  and  the  velocity  profiles 
reported  by  Burggraf  (28).  Furthermore,  the  difference  between  our  computed 
volume  flow  rate  and  that  predicted  by  the  empirical  characteristic  equation  of 
Denson  and  Hwang  (5)  for  a  single  screw  extruder  was  less  than  5%  for  a  wide 
range  of  pressure  rise  values  (Fig.  1b).  In  Fig.  1c,  the  temperature  distribution 
along  the  center-line  of  the  channel  for  a  temperature  dependent  shear 
viscosity  fluid  as  predicted  numerically  by  Gopalakrishna  et  al.  (29)  for  a 
parallel-plate  geometry  without  the  transverse  convection  terms,  is  compared 
with  our  computed  results,  also  ignoring  the  transverse  convection  terms.  It  is 
observed  that  the  agreement  is  indeed  very  good,  the  slight  difference  between 
the  results  can  be  attributed  to  the  fact  that  our  results  are  based  on  the  actual 
geometry  of  the  single  screw  extruder  with  the  curvature  included  to  give  rise  to 
slightly  higher  velocity  gradients  and  consequently  greater  viscous  dissipation 
effects  than  the  simplified  parallel  plate  based  calculations  of  Gopalakrishna  et 
al.  (29).  Other  essential  parts  of  the  code  have  been  similarly  compared  with 
numerical  results  available  in  the  literature  to  arrive  at  an  optimum  element  size 
distribution  for  accurate  numerical  predictions.  However,  experimental  results 
are  necessary  to  validate  the  model. 

The  only  experimental  study  which  was  aimed  at  determining  the 
temperature  distribution  in  the  extruder  was  carried  out  by  Esseghir  and  Sernas 
(30,31),  who  obtained  experimental  results  of  temperature  distribution  along  the 
center-line  of  the  channel  of  a  single  screw  extruder.  The  double  flighted  screw 


element  used  had  a  diameter  of  30.6mm,  channel  depth  of  0.31  Rs,  depth  to 
width  ratio  of  0.40,  a  helix  angle  of  16.2°  and  a  pitch  of  28mm.  There  were  three 
heat  transfer  sections  along  the  axis  of  the  barrel  with  the  barrel  temperature  in 
each  section  maintained  at  a  constant  value.  The  fluid  used  was  a  corn  syrup 
with  a  temperature  dependent  shear  viscosity,  the  rheological  and 
thermophysical  properties  of  which  were  determined  (30,31). 

The  first  set  of  experimental  temperature  measurements  chosen  for 
comparison  (Fig.  2)  were  obtained  with  the  first  barrel  section  maintained  at 
22oC,  the  same  temperature  as  that  of  the  fluid  at  the  inlet  of  the  extruder.  The 
solid/fluid  interface  temperatures  at  the  second  and  third  sections  were  kept  at 
11.90c.  In  the  second  set  of  temperature  measurements,  results  of  which  are 
shown  in  Fig.  3a,  the  barrel  wall  temperatures  for  all  the  three  heat  transfer 
sections  were  maintained  at  12.2oC,  while  the  inlet  temperature  of  the  fluid  was 
20.30c.  On  the  basis  of  the  geometrical  parameters  of  the  screw  element,  a 
finite  element  mesh  was  constructed.  The  required  dimensionless  parameters 
in  the  equations  of  conservation  of  mass,  momentum,  and  energy  were 
obtained  from  the  operating  conditions,  thermophysical  and  rheological 
properties.  For  both  experimental  runs,  the  value  of  PeuxH^/W  was  determined 

to  be  approximately  120,  which  suggests  that  transverse  convection  term  may 
not  be  negligible.  Rgures  2  and  3a  show  the  comparison  between  the  predicted 
and  experimental  temperature  distributions  with  the  origin  as  the  screw  root 
when  the  transverse  convection  terms  are  included  in  the  energy  equation.  The 
agreement  is  excellent  and  the  differences  are  within  typical  experimental  and 
numerical  error  ranges.  Furthermore,  the  temperature  distributions  confirm  that 
the  adiabatic  condition  at  the  screw  root  is  appropriate. 

If  the  transverse  convection  terms  in  the  equation  of  energy  were 
neglected.  Figs.  2  and  3a  also  indicate  that  the  numerical  predictions  would  be 
grossly  in  error.  Therefore,  this  simple  experiment  demonstrates  that  the 
commonly  employed  assumption  of  negligible  transverse  convection  in  single 
screw  extrusion  may  be  inappropriate.  Of  course,. for  a  sufficiently  long  extruder 
and,  negligible  viscous  dissipation,  as  we  approach  the  exit,  all  the  temperature 
values  will  approach  that  at  the  wall.  This  explains  why  the  difference  in 
temperature  values  between  the  predictions  with  and  without  transverse 
convection  terms  in  the  second  experiment  at  z*  of  48.  is  less  than  that  for  the 


first  6xp0rim0nt  at  z  =20.  Wh0n  th0  nnod0l  pr0dictions  ar0  compar0d  for  th0 
S0cond  0xp0rim0nt  at  z  of  20.,  Fig.  3b  shows  a  groator  difforanca  batwaan  tha 
predictad  temparature  valuas  with  and  without  transvarsa  convaction  tarms.  On 
tha  othar  hand,  the  corresponding  experimental  values  of  the  total  pressure 
generated  in  the  extruder  for  the  two  experimental  runs  were  21.0  and  20.4  bars 
while  the  numerical  predictions  are  15.4  and  16.3  bars  respectively.  This 
agreement  should  be  considered  to  be  satisfactory  noting  that  the  experimental 
operating  conditions  and,  the  geometry  of  the  extruder  were  such  that  the 
volume  flow  rate  versus  pressure  gradient  relationship  exhibits  a  steep  slope.  A 
5%  experimental  error  in  volume  flow  rate  would  account  for  the  differences 
obtained  and,  errors  of  this  magnitude  in  volume  flow  rate  measurements  are 
not  unexpected  considering  the  method  of  acquisition  of  the  data  (31). 

Since  the  Peclet  number  (3500)  for  the  experimentai  conditions  used 
(30,31)  is  not  high  enough,  the  use  of  upwinding  was  not  necessary.  However, 
for  high  Peclet  numbers,  the  use  of  upwinding  may  be  necessary.  For  example, 
with  the  Peclet  number  at  28,000,  l=igs.  4a-b  show  that  the  temperature 
distributions  obtained  with  and  without  the  use  of  SUPG  method  are  totally 
different.  This  value  of  Pe  is  based  on  the  operating  conditions  and  physical 
properties  of  Table  1,  which  are  in  the  range  we  typically  utilize  in  our  own 
experiments  on  industrial  scale  extruders.  The  unwound  screw  channel  profile 
is  that  of  a  screw  element  (Baker-Perkins  50.8mm  machine  )  with  a  channel 
depth  of  0.44Rs  and  depth  to  width  ratio  of  0.5.  It  is  obvious  that  at  this  Peclet 
number,  the  numerical  predictions  without  SUPG  are  inaccurate  since  the 
dimensionless  temperature  not  only  oscillates  but  also  violates  the  physics  of 
the  problem  which  precludes  a  temperature  value  less  than  that  at  the  inlet  (0o 
=  0.0).  Such  negative  values  of  temperature  have  also  been  reported 
elsewhere,  (32),  when  the  central  difference  scheme  (the  equivalent  of  Bubnov- 
Galer1<in  method)  is  implemented  in  the  discretization  of  the  convection  terms  of 
the  energy  equation.  When  the  Peclet  number  is  increased  to  2.8x1  O^,  the 
temperature  oscillations  become  more  pronounced  without  the  use  of  SUPG  as 
evidenced  from  Figs.  5a-b.  All  the  results  to  be  presented  in  the  next  section 
pertain  to  the  conditions  in  Table  1  and  other  necessary  parameters  are 
specified  in  the  figures. 

The  model  can  be  used  to  assess  the  importance  of  wall  slip  on 
development  of  temperature  distributions  in  extrusion  flows.  The  effects  of  the 


slip  parameter  on  the  dimensionless  temperature  distribution  are  depicted  in 
Figs.  6a-b.  At  the  down  channel  location  of  z*  =  20.0,  the  value  of  the  minimum 
dimensionless  temperature  within  the  channel  decreases  as  the  slip 
parameter,  p  ,  increases.  Wall  slip  essentially  reduces  the  local  temperature. 
The  corresponding  channel  direction  velocity  distributions  are  shown  in  Figs. 
7a-b.  For  the  same  volume  flow  rate,  as  the  slip  parameter  increases,  the  size 
and  strength  of  the  backflow  region,  and  the  maximum  downchannel  fluid 
velocity  decrease. 

The  development  of  the  dimensionless  pressure  gradient  along  the 
channel,  presented  in  Fig.  8,  shows  that  at  any  channel  location,  the  pressure 
gradient  generated  decreases  as  the  slip  parameter  Increases  for  a  fixed 
volume  flow  rate.  Wall  slip  reduces  the  velocity  gradients  thereby  reducing  the 
pressure  gradient.  In  addition,  very  close  to  the  entrance  region,  the  high 
velocity  gradients  result  in  an  increase  in  the  pressure  gradient.  Subsequently, 
a  decrease  in  viscosity  arising  from  the  temperature  rise  accounts  for  the 
decrease  of  the  pressure  gradient  down  the  channel.The  effect  of  the  velocity 
gradient  reduction  on  the  temperature  development  will  be,  as  expected,  to 
decrease  the  viscous  dissipation  effects  and  hence  result  in  lower 
dimensionless  bulk  temperature.  This  is  confirmed  by  the  results  in  Fig.  9  which 
show  a  decreased  value  of  dimensionless  bulk  temperature,  6b,  as  the  slip 
parameter  increases,  which  was  also  suggested  by  the  temperature  profiles 
reported  in  Figs.  6a-b. 

CONCLUSION 

A  non-isothermal  model  of  the  single  screw  extrusion  processing  of . 
generalized  Newtonian  fluids  which  includes  the  transverse  convection  terms  in 
the  energy  equation  has  been  developed  using  the  Streamline  Upwind/Petrov- 
Galerkin  (SUPG)  method.  The  numerical  predictions  of  the  code  are  In  excellent 
agreement  with  experimental  temperature  measurements  and  other  numerical 
results  available  In  the  literature.  The  capabilities  of  the  model  were  also 
demonstrated  by  investigating  the  effect  of  the  wall  slip  boundary  condition, 
commonly  observed  with  a  variety  of  materials  including  gels  and  concentrated 
suspensions.  This  model  should  provide  a  better  predictive  understanding  of 
the  single  screw  extrusion  process. 
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Table  1 


OPERATING  CONDITIONS  AND  PROPERTIES _ 

V  (V(moO)^))  =  0.05  (dimensionless)  Radius  of  barrel  (R^)  =  0.025m 

Hb*  (Hbco)  =  80  (dimensionless) _ Helix  angle  0  =  17.61° _ 

to  =  4.2  3-1 _ ___p=J400_kgAn^_ _ 

Ty=  6,650.  Pa _ C=2000J/kgK 

mo  =  35,000.  Pa-s  _ k  =  0.27W/mK _ 

n  =  1-0 _ _ _ |c‘(Tb-To)=2.0 _ 

Pe  =  28,000. 


List  of  Figures 


Figure  1:  Comparison  of  model  predictions  of  present  work  with  the  numerical 
results  of ;  (a)  Burggraf  (28),  (b)  Denson  and  Hwang  (5),  and  (c)  Gopalakrishna 
et  al.  (29). 

Rgure  2:  Comparison  of  numerical  results  of  temperature  distribution,  with  and 
without  transverse  convection  terms,  with  experimental  temperature  distribution 
at  z*=20.0  available  from  (30,31). 

Figure  3:  Comparison  of  numerical  results  of  temperature  distribution:  (a)  with  and 
without  transverse  convection  terms,  with  experimental  temperature  distribution 
3t  2*=48.0  available  from  (30,31),  (b)  with  and  without  transverse  convection 
terms  at  z*=20.0. 

Rgure  4:  Comparison  of  dimensionless  temperature  distributions  with  and  without 
SUPG  at  a  dimensionless  downchannel  location,  z*=10.0,  Pe=28,000,  and 
p*=0.0;  (a)  with  SUPG,  (b)  without  SUPG. 

Rgure  5:  Comparison  of  dimensionless  temperature  distributions  with  and  without 
SUPG  at  a  dimensionless  down  channel  location,  z*=20.0,  Pe=2.8x105,  and 
p*=0.1 ;  (a)  with  SUPG,  (b)  without  SUPG. 

Rgure  6:  Effect  of  slip  pararneter  on  the  dimensionless  temperature  distribution  0  at 
z*=20.0  for  Q*=.053;  (a)  p*=0.0,  (b)  p*=0.1 0. 

Figure  7:  Effect  of  slip  parameter  on  the  dimensionless  downchannel  velocity 
distribution  u^  at  z*=20.0  for  QV053;  (a)  p*=0.0,  (b)  p*=0.10. 


Figure  8:  Effect  of  slip  parameter  on  the  dimensionless  pressure  gradient  dpmVdz* 
forQ*=.093. 


Rgure  9:  Effect  of  slip  parameter  on  the  dimensionless  bulk  temperature  Bh  for 
QV093. 
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DEGREE  OF  MIXING  ANALYSES  OF  CONCENTRATED 
SUSPENSIONS  BY  ELECTRON  PROBE 
AND  X-RAY  DIFFRACTION 

R.  Yazici  and  D.  M.  Kalyon 

Highly  Filled  Materiai-s  Institute,  Stevens  Institute  or  Technolooy,  Hoiioken,  New  .Ieiisev  070:!0 

INTRODUCTION 

One  of  the  most  challenging  aspects  of  any  mixing  operation,  where  two  or  more  iden¬ 
tifiable  components  are  brought  together,  Is  the  characterization  of  the  state  of  the  mixture 
i.e.,  the  degree  of  mixedness  or  the  “goodness”  of  mixing.  In  the  nondilTusive  mixing  of  a 
viscous  polymeric  binder  with  solid  components,  the  complete  description  of  the  state  of 
the  mixture  would  require  the  specification  of  the  sizes,  shapes,  orientations  and  the  positions 
of  the  ultimate  particles  of  the  components. 

The  techniques  applied  to  characterize  the  degree  of  the  mixedness  of  a  composite  would 
involve  either  direct  or  indirect  methods.  Indirect  methods  include  the  determination  of 
some  representative  property  including  transmissive,  reflective,  resistivity,  ultrasound  and 
rheological  properties.  Such  techniques  are  ideal  especially  for  product  quality  control,  where 
the  relationship  between  the  measured  property  and  the  ultimate  property  of  interest  is 
apriori  known.  However,  although  such  techniques  are  useful,  they  offer  little  insight  into 
the  mechanisms  and  efficiency  of  the  mixing  operation. 

On  the  other  hand,  direct  measurement  techniques  can  provide  detailed  information  on 
the  process  and  the  development  of  the  microstructure  of  the  composite.  The  simplest  tech¬ 
nique  which  can  be  utilized  to  determine  the  mechanisms  of  mixing  in  a  batch  or  continuous 
mixer  is  the  injection  of  a  dye  into  a  transparent  fluid  and  then  to  follow  the  distribution 
and  the  interfacial  area  growth  of  the  nondiffusive  tracer.  Such  techniques  were  used  by 
Bigio  and  Erwin'  and  Bigio  et  al.^  to  analyze  the  mixing  mechanisms  occurring  in  tangential 
and  partially-full  co-rotating  twin  screw  extruders.  Ottino'  has  employed  the  flow  visual¬ 
ization  technique  to  investigate  the  mechanisms  of  mixing  in  various  time  dependent  flows, 
which  may  be  governed  by  chaotic  dynamics.  On  the  other  hand,  Morikawa  et  al*  have 
followed  the  spatial  distribution  of  color  incorporated  elastomers  to  gain  insight  into  the 
mixing  occurring  in  internal  mixers.  The  dispersive  mixing  aspects,  in  conjunction  with 
simple  shear  flows,  where  the  imposed  stresses  can  be  monitored,  were  also  characterized 
employing  the  cone  and  plate  flow  of  a  transparent  oil  containing  carbon  black  agglomerates 
(Rwei  et  al.^). 

If  the  materials  of  interest  are  opaque  or  if  transparent  barrel  sections  cannot  be  built, 
the  mechanisms  of  mixing  can  be  studied  through  "post-mortem”  analysis.  In  this  technique, 
generally  a  distinguishable  tracer  is  added  into  a  mixer  in  a  step  or  pulse  fashion.  After  a 
certain  duration  of  mixing,  the  mixture  is  systematically  removed  and  sectioned  into  parts 
to  allow  the  investigation  of  distributive  and  dispersive  mixing  aspects.  David  and  Tadmor," 
have  used  this  technique  for  probing  the  mechanisms  of  laminar  mixing  occurring  in  co¬ 
rotating  disk  processors  and  to  follow  dispersive  mixing  of  rubber  and  carbon  black  for¬ 
mulations.  Gokbora^  has  utilized  pigmented  polyethylene  granules  to  investigate  the  state 
of  mixedness  occurring  in  single  screw  extrusion  processing.  Kubota  et  al."  have  investigated 
the  flow  mechanisms  and  mixing  of  rubber  compounds  in  single  screw  extruders  based  on 
the  post-mortem  analysis  of  tracer  incorporated  rubber  strips  upon  screw  pulling.  Kalyon 
and  co-workers”"'  have  employed  color  incorporated  thermoplastic  elastomers,  followed  by 
computerized  image  analysis,  to  investigate  the  distributive  mixing  of  thermoplastic  elas¬ 
tomers  in  the  regular  flighted  and  kneading  disc  elements  of  twin  screw  extruders. 
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On  the  other  hand,  the  rapid  advent  of  imaging  and  sensing  technology  has  facilitated 
the  introduction  of  various  powerful  techniques,  including  tlie  magnetic  resonance  imaging 
and  x-ray  based  techniques,  to  the  analyses  of  opaque  mixtures.  Kalyon  et  al. ' '  have  employed 
magnetic  resonance  imaging,  wide-angle  x-ray  diffraction  and  x-ray  radioscopy  to  characterize 
the  amount  of  air  entrained  into  composite  suspensions  with  elastomeric  binders  during 
extrusion  processing,  which  is  related  principally  to  the  degree  of  1111  profile  In  the  extruder. 
Techniques  such  as  x-ray  computed  tomography  were  applied  to  the  study  and  the  quality 
control  of  advanced  composites, rocket  motors,'’  precision  castings,”'"’’  ceramics'"  and 
assembled  structures'^  but  have  not  been  applied  to  the  study  of  the  state  of  mixedness. 
The  applications  of  various  nondestructive  evaluation  techniques  in  the  tire  industry  are 
reviewed  by  Trivisonno.'"  It  may  be  envisioned  that  such  powerful  techniques,  including 
magnetic  resonance  imaging  and  computed  tomography,  with  three  dimensional  and  real¬ 
time  analysis  capabilities,  will  be  utilized  widely  for  on-line  process  and  product  control  of 
various  mixing  and  processing  flows  in  tlie  near  future  in  many  industries. 

In  this  study  the  energy  dispersive  x-ray  and  wide  angle  x-ray  diffraction  techniques  are 
being  introduced  into  the  field  of  goodness  of  mixing  measurements.  The  applicability  of 
these  two  techniques  to  the  quantitative  characterization  of  the  distributive  mixing  occurring 
in  continuous  mixers  and  the  scale  of  examination  effects  has  been  demonstrated  with  a 
highly  filled  suspension  involving  an  elastomeric  binder,  i.e.,  hydroxyl  terminated  polybu¬ 
tadiene,  mixed  with  two  different  solid  fillers. 

EXPERIMENTAL  PROCEDURES 
MATERIAL  PROCESSING  AND  SAMPLE  PREPARATION 

In  this  study,  the  concentrated  suspension  consisted  of  a  polymeric  matrix  and  76.5 
volume  percent  of  solids.  The  matrix  was  hydroxyl  terminated  polybutadiene  (HTPB)  and 
an  antioxidant.  The  fillers  were  aluminum  (Al)  and  ammonium  sulfate  (AS)  in  powder  form. 
Mean  particle  size  of  aluminum  was  20  microns.  Ammonium  sulfate  was  bimodal  in  particle 
size  with  75%  coarse  particles  and  25%  fine  particles.  The  mean  particle  sizes  of  the  coarse 
and  fine  fractions  of  ammonium  sulfate  were  200  and  20  microns,  respectively. 

An  APV  MP50  twin  screw  extruder  with  a  screw  length-to-diameter  ratio  of  15  was 
employed  for  processing  of  the  suspension  samples.  The  extruder  was  equipped  with  a  Zenith 
gear  pump  for  feeding  of  the  liquid  ingredients,  a  K-Tron  volumetric  feeder  for  the  feeding 
of  aluminum  and  an  Acrison  loss-in-weiglit  feeder  for  the  feeding  of  ammonium  sulfate.  The 
moisture  content  of  ammonium  sulfate  was  reduced  to  below  0.r^>,  before  it  was  placed  into 
the  hopper  .of  the  loss-in-weight  feeder.  This  feeder  was  kept  sealed  at  all  times.  The  con¬ 
ditions  of  the  suspension  preparation,  including  screw  configuration  and  screw  speed  (43 
rpm)  were  kept  the  same  in  all  of  the  experiments.  A  vacuum  port  was  installed  and  employed 
in  conjunction  with  specially  designed  devolatilization  elements,  “ball  wings”,  to  maximize 
the  surface-to-volume  ratio  of  the  suspension,  for  the  satisfactory  removal  of  the  air  entrained 
at  the  earlier  stages  of  processing.  The  devolatilization  section  was  sealed  by  the  die  in  the 
upstream  direction  and  a  series  of  reversely  configured  fully-flighted  conveying  screw  ele¬ 
ments  in  the  downstream  direction.  A  Welch  Duo-Seal  vacuum  pump  was  used  to  remove 
the  air  from  the  suspension  through  a  vacuum  port. 

The  samples  were  collected  from  the  exit  of  the  die  under  steady-state  conditions.  Sample 
cups  with  dimensions  of  45  X  15  mm’  and  with  1  mm  depth  were  filled  flush,  with  the 
suspension.  The  samples  were  kept  frozen  at  all  times  to  prevent  segregation  effects,  except 
during  the  duration  of  the  microstructural  analyses. 

CHARACTERIZATION  OF  MIXING  INDICt-S 

Various  mixing  indices,  which  characterize  the  goodness  of  mixing  of  the  suspension 
samples,  were  determined  on  the  basis  of  the  relative  volume-fraction  measurements  of  the 
two  filler  materials  by  tliree  principal  analytical  methods:  scanning  electron  microscopy 


/■n  n’-n-qt  07:56: is  i.P:  rcpT  —  noo? 


MIXING  ANALYSIS 


0039 


pEM),  electron-probe  energy-dispersive  x-ray  analysis  (EDX)  and  wide-angle  x-ray  dilTrac- 
tometry  (XRD).  Mixing  indices  were  determined  at  different  scales  of  examination  by  varying 
the  sample  area  through  alteration  of  the  size  of  the  electron  and  x-ray  probes  by  several 
orders  of  magnitude.  The  ratio  of  the  volume  fractions  of  the  two  fillers,  aluminum  (Al)  and 
ammonium  sulfate  (AS),  itiA,/4>As,  was  utilized  as  a  measure  of  degree  mixedness  of  the  sample. 

In  general,  the  quantitative  description  of  the  mixing  quality  or  goodness  of  mixing  of  a 
given  mixture  can  be  developed  by  comparison  of  the  state  of  the  mixture  to  the  most  complete 
mixing  state  attainable.  "  This  complete  mixing  corresponds  to  statistical  randomness  of 
the  ultimate  particles  of  the  ingredients  being  mixed. “ 

If  one  makes  N  measurements  of  concentration  c.  of  one  of  the  components,  then  the 
mean  concentration  may  be  calculated  according  to 


-  ^  V 


(1) 


where  c  should  not  differ  significantly  from  the  overall  concentration  of  the  component 
unless  a  faulty  sampling  technique  is  used.  The  difference  between  c  and  ^  decreases  as  the’ 
finite  number  N,  of  the  characterized  samples,  is  increased.  The  measured  concentration 
values  of  the  minor  component  also  depend  on  the  sample  size.  These  values  approach  the 
overall  concentration  of  the  minor  component  <l>  as  the  sample  size  is  increased. 

measure  of  the  homogeneity  of  a  mixture  is  the  extent  to  which  the  concentration 

^bpTnri  fom  the  mean  concentration. 

lance  s  ,  arising  from  the  individual  concentration  c,  measurements,  provides  such 
an  index  to  quantitatively  assess  the  degree  of  mixedness.  The  variance  is  given  by 


s'  = 


•  2  (C(  -  cf. 


(2) 


are  clniT'M I"’'!!*"'’  concentration  c.  values  which 

variance  o”cmiJs"irM  a  homogeneous  system.  The  maximum 

^ariance  occurs  if  the  components  are  completely  segregated.  Maximum  variance  is  given 


sii  =  c(  1  -  C  ). 


(■i) 


If  the  variance  is  normalized  to  its  maximum  value,  the  resulting  parameter  is  called  the 
intensity  of  segregation,  /,„.  This  is  given  by 


/  =5!  =  _iL_ 

So  c(  I  -  C  ) 


(4) 

important  index  in  evaluating  goodness  of  mixing  or 
distribution.  The  intensity  of  segregation  ranges  from  unity,  for  completely  segregated  sys 
tern,  to  zero,  for  a  homogeneous  system.  segregacea  sys 

In  this  study,  the  mean,  variance  and  intensity  of  segregation  of  the  relative  concentration 
fractions  of  the  two  fillers  (aluminum  and  ammonium  sulfate),  <l>A,/<i>A’<,  were  determined  to 
assess  the  goodness  of  mixing.  The  overall  value  of  the  relative  volult ncentratS^^^^ 
of  aluminum  and  ammonium  sulfate,  was  0.20  ±  0.02,  as  specified  during  processing. 

ANALYTICAL  TECHNIQUES 

Electron  microscopy  coupled  with  energy-dispersive  x-ray  analysis  is  a  powerful  tool  in 

nulti  Z^  f corresponding  elemental  distributions  in 
u Iti-phase  materia  s.  In  this  study  the  two  techniques  were  utilized  in  synergy  to  obtain 
b  )th  qualitative  and  quantitative  information  on  the  state  of  mixedne.ss.  First  the  spatial 

order  to  distinguish  the  aluminum  and  ammonium  sulfate  particles  from  each  other,  el- 
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emental  mappings  of  aluminum  and  sulfur  were  carried  out  with  energy-dispersive  x-ray 
analysis  at  the  same  locations,  where  the  secondary-electron  images  were  taken.  The  identity 
of  each  particle  in  the  SEM  images  was  thus  determined  for  further  evaluation. 

Chemical  analysis  of  materials  in  the  scanning  electron  microscope  is  performed  by  mea¬ 
suring  the  energy  and  intensity  distribution  of  the  x-ray  signal  generated  by  a  focused 
electron  beam.  A  sufllciently  energetic  electron  beam  may  eject  an  inner-shell  electron,  leaving 
the  atom  in  an  ionized  or  excited  state  with  a  vacancy  in  one  of  the  shells.  During  subsequent 
deexcitation  of  the  atom,  electron  transition  occurs  from  an  outer  shell  to  fill  this  vacancy. 
The  transition  involves  a  change  in  energy.  The  excess  energy  which  the  excited  atom  contains 
can  be  released  in  the  form  of  a  photon  in  the  x-ray  range  of  the  electromagnetic  spectrum. 
Since  the  energy  of  the  emitted  x-ray  is  related  to  the  difference  in  energy  between  the 
sharply  defined  levels  of  the  particular  atom,  it  is  referred  to  as  a  characteristic  x-ray.  The 
specific  energies  of  the  characteristic  x-ray  peaks  for  each  element  has  been  tabulated  and 
can  be  used  to  evaluate  the  unknown  materials. 

The  energy-dispersive  x-ray  analyzer,  shown  schematically  in  Figure  1,  employs  a  solid- 
state  detector  to  acquire  tlie  total  spectrum  of  characteristic  x-rays  from  0.75  to  20  keV  (or 
more).  Large  numbers  of  x-ray  photons  can  be  counted  at  all  energy  levels  with  the  aid  of 
a  multi-channel  analyzer.  This  allows  for  the  rapid  evaluation  of  numerous  elements  in  the 
sample. 

The  intensity  of  the  characteristic  x-radiation  is  proportional  to  the  concentration  of  the 
particular  element  in  the  sampling  volume.  Intensity  values,  therefore,  can  be  utilized  for 
quantitative  analysis.  The  number  of  characteristic  x-ray  photons  per  second  (A^()  detected 
from  an  element  is  given  by: 

N,=Jk,<l,r,  (5) 

where  J  is  tl»e  incident  electron  flux,  4i,  is  the  concentration  of  the  element  t,  k,  is  a  propor¬ 
tionality  factor  that  includes  ionization  cross-section,  fluorescent  yield,  and  detection  effi¬ 
ciency. 
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With  computerized  modern  instruments  the  values  are  readily  calculated  and  quanti¬ 
tative  analysis,  accurate  to  0.5  weight  percent,  is  achieved  for  elements  with  an  atomic 
number  higher  than  sodium.^'  The  absolute  determination  of  the  elemental  concentration, 
however,  requires  calibration  with  standard  samples,  especially  when  light  elements  are 
present  in  the  sample.  Determination  of  relative  concentrations  of  two  elements  in  a  specimen 
can  be  reduced  to: 

Nj 

where  fcn  is  the  relative  proportionality  factor  for  the  two  elements  in  a  given  mixture. 

in  this  study  the  relative  atomic  fractions  of  the  two  elements  aluminum  and  sulfur  were 
measured  by  energy-dispersive  x-ray  analysis.  Variation  of  the  relative  atomic  fractions  of 
the  two  elements  was  utilized  to  assess  the  degree  of  distributive  mixedness.  In  order  to 
calibrate  the  measured  values  standard  samples  were  prepared  by  mixing  controlled  amounts 
of  aluminum  and  ammonium  sulfate  powders.  According  to  this  calibration  the  relative 
atomic  fraction  proportionality  factor  /croi.  was  1.21,  i.e., 

N^,/Ns=  (7) 

where  <t>M  and  <t>f  refer  to  the  atomic  concentrations  of  aluminum  and  sulfur,  respectively. 

The  relative  weight  fraction  and  relative  volume  fraction  values  for  aluminum  and  am¬ 
monium  sulfate  components  were  then  calculated  from  the  relative  atomic  fraction  of  alu¬ 
minum  and  sulfur.  The  conversion  factor  from  relative  atomic  to  relative  weight  fraction 
was  4.898  and  from  relative  atomic  to  relative  volume  fraction  was  7.47.  Molecular  and 
crystal  structure  parameters'^  of  the  two  components  were  used  in  determining  these  con¬ 
version  factors.  Finally,  the  relative  volume  fraction  proportionality  factor  ,  was  9.04, 
i.e., 

Na,/Ns’=9M4>7,‘/<1>"^-  (8) 

Energy-dispersive  x-ray  measurements  were  carried-out  with  a  series  of  successively 
smaller  probes  at  sites  chosen  systematically,  based  on  a  grid  system.  A  JOEL  840  (40  A 
resolution)  scanning  electron  microscope  with  Tractor  Northern  TN  5500  energy-dispersive 
x-ray  analysis  system  was  used  in  these  studies.  Incident  beam  voltage  was  set  at  15  kV. 
The  take-off  angle  of  the  x-rays  (the  angle  of  the  x-ray  detector  with  respect  to  the  sample 
surface)  was  set  at  40  degrees.  The  window  thickness  was  7 .5  microns.  Multiple  probe  sizes 
were  used  in  these  experiments  (80  mm^,  9  mm*,  1  mm*,  0.1  mm*  and  0.01  mm*)  to  determine 
the  distributive  mixing  index,  i.e.,  intensity  of  segregation,  as  a  function  of  the  scale  of 

examination.  _  . 

X-ray  diffractometry  has  been  successfully  applied  for  both  qualitative  and  quantitative 
phase  analysis  in  multi-phase  materials.*^’**  Upon  irradiation  with  x-rays,  a  given  substance 
produces  a  characteristic  diffraction  pattern.  Diffraction  is  essentially  a  scattering  phenom¬ 
enon  in  which  a  large  number  of  atoms  in  the  crystalline  material,  arranged  periodically  on 
a  lattice,  scatter  the  x-rays  in  phase.  These  phase  relations  are  such  that  destructive  inter¬ 
ference  occurs  in  most  directions  of  scattering;  but  in  a  few  directions  constructive  inter¬ 
ference  takes  place  and  diffracted  beams  are  formed.  According  to  Bragg’s  Law,  the  diffracted 
beams  make  a  28  angle  with  respect  to  the  Incident  beam,  i.e. 

\  =  2d  Sin  9,  (9) 

where  X  is  the  wavelength  of  the  x-rays  of  the  incident  beam  and  d  is  the  spacing  between 
the  atomic  (lattice)  planes  of  the  crystalline  material.  Since  the  d-spacings  between  atomic 
planes  and  their  distribution  in  space  (the  crystal  structure)  is  unique  for  each  material, 
the  angular  distribution  of  the  diffraction  peaks  and  their  intensities  (the  diffraction  pattern) 
is  also  unique  for  a  particular  material.  Diffraction  patterns  are  obtained  with  diffractometers 
equipped  with  suitable  detectors.  Qualitative  analysis  by  x-ray  diffraction  is  accomplished 
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by  identification  of  the  particular  diffraction  pattern  of  a  substance  from  the  standard  dif¬ 
fraction  tables. 

The  particular  advantage  of  x-ray  diffraction  analysis  is  that  it  discloses  the  presence 
of  a  substance,  as  that  substance  actually  exists  in  the  sample,  and  not  In  terms  of  its  con¬ 
stituent  chemical  elements.  If  the  sample  contains  more  than  one  compound  or  phase  that 
constitute  the  same  chemical  elements,  all  these  compounds  are  disclosed  by  diffraction 
analysis.  Quantitative  analysis  is  also  possible,  because  the  intensity  of  the  diffraction  pattern 
of  a  particular  phase — in  a  mixture  of  phases — depends  on  the  concentration  of  that  phase 
in  the  mixture.  The  relation  between  the  Integrated  intensity  4  and  the  volume  fraction  4>r 
of  a  phase  is  generally  nonlinear  since  diffracted  intensity  depends  strongly  on  the  absorption 
coefficient  of  the  mixture,  Mm,  vvhich  itself  depends  on  the  concentration.  For  a  two-phase 
material,  (with  absorption  coefficients  /ii  and  m2  for  the  individual  phases)  the  absorption 
coefficient  for  the  mixture  becomes: 


Mm  =  •  (10) 

Tlie  integrated  intensity  from  phase  1  is  then  given  by''*: 

I\  —  (11) 

where  K,  is  a  constant  that  depends  on  the  material  and  the  incident  beam  used  but  not  on 
the  concentration.  The  ratio  of  intensities  from  phases  1  and  2,  however,  is  independent  of 
Mm  and  varies  linearly  with  concentration: 

hUi  =  (12) 

In  this  study  the  ratio  of  relative  volume  fractions  of  aluminum  and  ammonium  sulfate, 
was  calculated  from  the  integrated  intensity  ratio  4i/4,v.  These  measurements  were 
carried  out  utilizing  particular  crystal-planes  of  the  two  filler  materials.  Relatively  high 
number  of  crystal-plane  reflections  were  evaluated  in  order  to  eliminate  texture  effects  and 
to  increase  accuracy.  Integrated  intensities  of  the  first  two  reflections  (ill  and  200)  of 
aluminum  and  the  first  five  reflections  (002,  Oil,  102,  111  and  200)  of  ammonium  sulfate 
were  utilized  for  calculations.  The  value  of  the  constant  was  taken  as  0.489  for  (111) 

planes.  This  value  was  calculated  according  to  the  i.,m  values  listed  by  JCPDS.^''  The 

contribution  of  the  absorption  of  tlie  matrix  to  the  absorption  coelTicient  of  the  mixture  was 
less  than  3%  and,  thus,  it  was  ignored.  The  relative  volume  fraction  <l>M/'i>As  values  were 
further  calibrated  with  measurements  carried-out  on  standard  samples. 

An  automated  GE  wide-angle  x-ray  diffractometer  was  used  for  these  studies.  Nickel 
filtered  Cu  K„  radiation  at  30  kV  and  15  mA  was  used.  The  scanning  speed  was  one  degree 
per  minute.  A  0.2  degree  receiving  slit  was  used  in  all  runs.  The  x-ray  probe  size  was  varied 
by  using  3  degree  and  0.4  degree  primary  beam  slits  and  a  8  mm  high  window.  The  x-ray 
probe  sizes  used  were  20  mm  X  8  mm  (160  mm^)  and  1  mm  X  8  mm  (8  mm"'),  respectively, 
at  a  Bragg  angle,  9,  of  20  degrees. 

RESULTS  AND  DISCUSSION 

Typical  SEM  photomicrographs  of  the  concentrated  suspension  sample  are  shown  in 
Figure  2  at  two  different  magnifications.  The  smaller  light-colored  particles  are  aluminum 
and  the  larger  dark-gray  particles  are  ammonium  sulfate.  The  matrix  is  the  dark  colored 
material  found  in  between  the  particles.  Some  regions  which  are  richer  in  either  aluminum 
or  ammonium  sulfate  can  be  discerned.  However,  it  is  not  possible  with  the  available  SEM 
technology  to  distinguish  and  quantitatively  characterize  the  concentration  of  each  solid 
component  at  a  given  location.  The  typical  characteristic  dimension  of  grouping/clustering 
of  alike  particles  is  400  to  600  microns.  Pockets  of  matrix-rich  regions  are  not  observed  upon 
microscopic  examination. 

It  is  difficult  to  prepare  samples  with  the  image  quality  as  indicated  in  Figure  2(a  and 
b).  As  sliown  in  Figure  3  (panel  a),  usually  the  matrix  covers  the  surface  particles  and  l)liirs 
the  images.  Removal  of  the  matrix  from  the  surface  without  disturbing  the  particle  distri- 
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Fio.  2.-Scanning  electron  micrographs  at  two  magnincatlons:  (a)  00  x,  (b)  200  X, 

bution  is  not  an  easy  task.  Keeping  the  sample  under  vacuum  for  an  extended  time  or  selective 
etching  of  the  matrix  works  with  limited  success.  In  most  cases  a  matrix  layer  covers  the 
surface  and  prevents  the  identification  of  the  location  of  the  particles.  In  order  to  identify 
the  aluminum  and  ammonium  sulfate  particles  and  to  carry-out  the  quantitative 
mixing  indices,  the  elemental  mappings  of  aluminum  and  sulfur  were  carned-out  using  energy 
dispersive  x-ray  analysis,  EDX,  at  the  same  locations  of  the  scanning  electron  micrographs. 
These  elemental  mappings  are  shown  in  Figure  3,  panel  b  and  c,  respectively.  As  shown  with 
connecting  arrows  in  Figure  3,  the  EDX  technique  can  clearly  distinguish  ammonium  sulfate 

and  aluminum  particles. 
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Fin.  4. — Variation  of  the  relative  volume  fraction  of  the  two  solid  components,  aliirniinim  and  ammonium  sniratr, 
at  dilTerent  locations  of  the  mixture  as  a  function  of  scale  of  cxamitialioti  of  enerKy-disporsive  x>ray  analysis. 


The  results  of  the  quantitative  energy-dispersive  x-ray  analysis  of  the  relative  volume 
fractions  of  aluminum  and  ammonium  sulfate,  4>Ai/<f>As,  are  given  in  Figure  4.  The  statistical 
analysis  of  the  results  are  listed  in  Table  I.  The  variance  of  the  relative  volume  fraction  0,„/ 
(tiAs,  depends  on  the  probe  size.  The  largest  probe  size,  80  mm"*  gives  ri.se  to  the  minimum 
variation  in  the  volumetric  ratio  as  expected.  With  the  moderate  probe  sizes,  1  mm^  to  10 
mm^  range,  the  variance  of  the  volume  ratio  remains  the  same.  At  the  relatively  smaller 
probe  sizes  of  0.1  mm^  and  0.01  mm^,  the  deviation  of  the  relative  volume  fraction  of  the 
two  solid  components  from  the  mean  increases  considerably  with  decreasing  probe  size.  The 
marked  change  in  the  variance  of  the  relative  volume  fraction  values  between  probe  sizes 
0.1  mm^  to  1  mm^  is  in  good  agreement  with  the  qualitative  microstructural  features  observed 
in  scanning  electron  microscopy.  Clustering  of  the  particles  in  300  to  600  micron  diameter 
regions  should  give  rise  to  the  observed  scatter  in  the  energy-dispersive  x-ray  analysis  results. 


Table  I 


Statistical  Parameters  of  the  Relative  Volume  Fractions,  'f.ii/'t’A'!,  for  Different  EDX  Probe  Size.s 


Relative  volume 
fraction, 
parameters 

EDX  probe  size 

80  mm'* 

9  mm* 

1  mm* 

0.1  mm* 

0.01  mm* 

mean 

0.216 

0.215 

0.217 

0.217 

0.216 

.standard  deviation 

0.007 

0.024 

0.025 

0.034 

0.060 

!)!)%  confidence  interval 

±0.006 

±0.014 

±0.013 

±0.018 

±0.030 

intensity  of  segregation 

0.0003 

0.0034 

0.0037 

0.0068 

0.0281 

phi 
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Probes  smaller  than  0.01  mm^  in  size  were  not  used,  since  they  are  in  the  same  size  range 
as  the  particle  diameter  of  ammonium  sulfate. 

The  99%  confidence  intervals  of  the  energy-dispersive  x-ray  analysis  data  determined 
according  to  Student’s-t-distribution  are  also  Included  in  Table  1.  The  mean  of  the  data  for 
all  probe  sizes  is  around  0.215-0.217.  These  mean  values  agree  with  the  input  volumetric 
ratio,  thus  suggesting  the  absence  of  systematic  errors.  On  the  other  hand,  the  standard 
deviation  values  increase  from  0.007  at  the  probe  size  of  80  mm®  to  0.07  at  a  probe  size  of 
0.01  mm®,  indicating  that  the  deviation  from  the  mean  increases  by  about  one  order  of 
magnitude  going  from  the  largest  to  smallest  probe.  The  confidence  intervals  again  reflect 
the  same  observation,  with  an  increasing  range  around  the  mean  as  the  probe  size  is  reduced. 
Obviously,  if  the  probe  size  is  reduced  further,  one  would  probe  only  the  completely  seg¬ 
regated  regions  consisting  of  either  aluminum  or  ammonium  sulfate  single  particles. 

Which  probe  size  (scale  of  examination)  should  be  used  for  the  characterization  of  the 
goodness  of  a  given  mixture?  The  scale  of  examination  should  be  selected  depending  on  the 
application  and  the  relevant  ultimate  properties  of  the  mixture.  For  some  applications  the 
sample  area  may  be  as  large  or  larger  than  the  maximum  probe  size  used  in  this  analysis, 
i.e.,  80  mm®.  However,  for  other  applications  it  may  be  necessary  to  achieve  a  certain  degree 
of  mixedness  at  a  smaller  scale  of  examination. 

The  results  of  the  quantitative  x-ray  dilTraction  analysis  of  relative  volume  fractions  of 
aluminum  and  ammonium  sulfate,  by  employing  the  relative  integrated  intensities,  are  given 
in  Figure  5.  The  statistical  analysis  of  these  results  are  presented  in  Table  II.  With  x-ray 
diffractometry,  the  mean  of  the  relative  volume  fraction  values  of  aluminum  In  the  solid 
phase  4>mI4>as,  was  determined  to  be  0.216.  Determination  of  the  relative  integrated  intensity 
//x.?,  values  at  the  smaller  probe  size  (8  mm®)  was  found  to  be  complicated  by  the  preferred 
crystallographic  orientation  exhibited  by  the  ammonium  sulfate  particles.  This  preferred 
orientation  was  especially  evident  with  002,  200,  or  202  planes.  As  seen  in  Figure  3,  the 


Fkj.  r, — Variation  of  llic  rniativc  voliimi'  fraction  of  the  two  soliil  componi-nt.-i.  aluminum  and  ammoninm  .snifalc. 
at  dilTc'rpnt  locations  of  tlif  mixture  a.s  a  function  of  scale  of  examination  of  x-ray  flilTractomctry. 
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Table  II 

Statistical  Pabameters  of  the  Relative  Volume 
Fractions,  for  Different  XRD  Probe  Sizes 


Relative  volume  XRD  probe  size 


fraction, 

parameters  160  mm*  8  mm* 


mean 

0.216 

0.216 

standard  deviation 

0.012 

0.029 

99%  confidence  interval 

±0.012 

±0.017 

intensity  of  segregation 

f).0009 

0.0050 

ammonium  sulfate  particles  exhibit  relatively  flat  and  angular  surfaces,  which  can  be  pref¬ 
erentially  oriented  during  processing.  Since  the  samples  are  taken  from  random  locations 
and  at  random  orientations  from  the  extruded  suspension,  a  variety  of  preferred  orientation 
characteristics  were  observed  during  x-ray  diffractometry. 

The  x-ray  diffractometry  analysis  also  indicates  that  the  various  values  of  deviation 
from  the  mean  increase  with  the  decreasing  scale  of  examination.  However,  there  are  dif¬ 
ferences  between  the  wide  angle  x-ray  diffraction  and  energy-dispersive  x-ray  analysis  data 
that  emanate  from  the  applicable  depth  of  penetration.  The  electron  beam  used  in  energy- 
dispersive  x-ray  analysis  is  strongly  absorbed  at  the  surface  thus  the  x-ray  photons  that 
make  it  to  the  detector  come  mostly  from  the  top  few  microns  depth.  This  depth  is  smaller 
than  the  particle  sizes  of  both  solid  components  of  the  formulation.  However,  in  the  wide- 
angle  x-ray  diffraction  measurements  96%  of  the  information  is  derived  from  the  top  50  to 
100  microns  of  the  suspension,  depending  on  the  filler  type  and  the  diffraction  angle  0.  The 
penetration  depth  can  be  estimated  from; 

a;  =  3  sin  •  (13) 

The  penetration  depth  obtained  with  the  x-ray  diffraction  technique  is  greater  than  the 
electron  probe  measurements  and,  thus,  x-ray  diffraction  results  are  more  representative  of 
the  bulk. 

The  intensity  of  segregation  index  values  were  determined  from  the  energy-dispersive 
x-ray  and  x-ray  diffractometry  data,  as  shown  in  Tables  I  and  11.  The  intensity  of  segregation 
values  increase  with  the  decreasing  scale  of  examination  in  both  techniques.  The  intensity 
values  obtained  with  energy-dispersive  x-ray  analysis  and  x-ray  diffraction  are  close,  in¬ 
dicating  that  both  techniques  are  capable  of  characterizing  the  degree  of  distributive  mixed¬ 
ness  of  the  solid  ingredients. 

It  should  be  noted  that  the  applicabilities  of  the  two  techniques  to  various  processes  and 
materials  are  different.  Wide-angle  x-ray  diffraction  may  also  be  suitable  for  on-line  analysis 
of  microstructural  distributions  and,  thus,  product  quality  control.  However,  the  scale  of 
examination  cannot  be  reduced  beyond  a  certain  threshold  with  x-ray  diffraction,  when 
preferred  orientation  effects  complicate  the  intensity  ratio  measurements.  On  the  other  hand, 
the  energy-dispersive  x-ray  analysis  is  post-mortem  in  nature  and  requires  a  lengthy  analysis 
and  application  of  vacuum  during  characterization. 

With  a  conventional  EDX  detector,  as  used  in  this  study,  the  elemental  mappings  are 
restricted  to  elements  with  atomic  numbers  above  that  of  stidium.  However,  with  other 
currently  commercially  available  EDX  units,  with  windowless  detectors,  elements  with  atomic 
numbers  greater  than  or  equal  to  five  can  also  be  mapped.  On  the  other  hand,  the  x-ray 
diffraction  technique  is  applicable  to  all  .semi-cry.stalline  ingredients,  regardle.ss  of  their 
atomic  mass. 
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CONCLUSIONS 

Two  x-ray  based  techniques  involving  energy -dispersive  analysis  and  dilTractometry  were 
introduced  to  the  analyses  of  the  degree  of  mixedness,  i.e.  the  ‘‘goodness  of  mixing”  of 
concentrated  suspensions.  A  hydroxyl  terminated  polybutadene  matrix  was  mixed  with 
aluminum  and  ammonium  sulfate.  In  the  analysis,  the  ratio  of  the  relative  volume  fractions 
of  the  two  solid  components  was  used  as  the  basis  of  the  analytical  evaluation.  Both  char¬ 
acterization  techniques  are  capable  of  determining  the  relative  volume  fraction  of  the  two 
solid  components  as  a  representative  measurement  of  the  distributive  mixing  efficiency  and 
both  are  sensitive  to  the  scale  of  examination.  The  introduced  techniques  should  be  useful 
in  the  better  definition  of  the  degree  of  mixedness  as  well  as  in  resolving  differences  in 
mixing  efficiencies  of  various  mixers  used  in  processing  of  concentrated  suspensions. 
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